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DESCRIPTION 

CHARAC TERISTIC V ALU E IDENTIFIC A TION METHOD^ 
AND APPARATUS 



Technical Field 

The present invention relates to a characteristic value 
identification method and an apparatus therefor, and in particular to a 
characteristic value identification method and an apparatus therefor of 
a part composing a product. 

Background Art 

There have been known two types of prior art characteristic value 
identification methods and apparatuses therefor in the following: 

One is disclosed in the Japanese Patent Application No.6-140081, 
No.8- 18946, No.9-287161, and the like. This prior art method adopts 
means for identifying a natural frequency, a natural mode, and a 
modal damping ratio based on measured data of a frequency response 
function (FRF) obtained by a vibration test of a tested object (part). 

The other is a method called a characteristic matrix identification 
method. This method adopts means for identifying three kinds of 
characteristic matrices, i.e. a mass matrix, a damping matrix, and a 
stiffness matrix, which are coefficients of an equation of motion 
formulated based on "equilibrium of forces" rule in a dynamics system, 
from time history measured data of an exciting force and a response 
obtained by the vibration test of the tested object (part). 

Generally, a machine has a mechanical power source generating a 
kinetic energy in some way, and is composed of a large number of parts 
such as a mechanism which transmits the mechanical power and a 
mechanism which takes advantage of the mechanical power for the 
work. Also, with the parts composing the machine, parts based on a 
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physical unit system such as electric, mechanical, and fluid systems 
are organically united or combined. 

The fact is that models which reproduce such a large number of 
parts with many different theoretical backgrounds are individually 
5 modeled for the identification according to the physical unit system in 
which the parts are included or the application thereof, with a mutual 
theoretical relationship being neglected regardless of high generality of 
the individual part and mechanism. For this reason, it has been 
disadvantageous that the above-mentioned two types of prior art based 
10 on the techniques of the mechanical system can be applied only to a 
M= specific phenomenon of mechanical system vibration. 
H Also, in the prior art, the model (part) has been identified only by 

one dimension "strength (flow quantity)" for two dimensional 
Si quantities "quantity (potential quantity)" and "strength" which 

m 

15 prescribe energy. Namely, since both rules which govern each of the 

L two kinds of quantities composing energy are not applied, continuity of 

Of velocity and acceleration has not been represented on the identified 

M* model, so that it was disadvantageously difficult to unite the parts 

o 

M= with each other. 

20 Furthermore, since the prior art only expresses a structure of a 

model, i.e. a positional relationship of the characteristics composing 
the model but not the functions of the characteristics, structural 
phenomena which the characteristics are compounded to reveal, e.g. a 
natural frequency, a natural mode shape, or a characteristic matrix of 

25 an entire structure can only be made clear by the identification, so that 
it has been impossible to directly clarify the characteristic values 
governing the function of an object. 

It is accordingly an object of the present invention to provide a 
characteristic value identification method and an apparatus therefor 

30 which can develop models across different physical unit systems to be 
integrated into an entire model, similar to a product in which 



2 



individual parts are combined. 



Disclosure of Invention 

A basic process for identifying an internal characteristic value of 
5 a functional model relating to a product and a part (hereinafter both 
simply referred to as a part) in the present invention is shown in Fig.l. 
This diagram illustrates an identification method taking advantage of 
the functional model being a modeling method including a steady 
characteristic value and a transient characteristic value. 
10 First of all, the functional model modeled in a process 1 originally 

q includes the steady characteristic value and the transient 
Jf characteristic value. Accordingly, the functional model identification of 
P the present invention according to claim 1, these internal 

m characteristic values are separately carried out by a steady 

T 15 identification process 2 and a transient identification process 3. 

Namely for a basic process of the identification, the transient 
identification process 3 is executed after the steady identification 
process 2, since the steady characteristic value is free from 
interference of transient state and the transient state interferes with 
20 the steady characteristic value. 

The steady identification process 2 further includes the following 
processes according to claims 2-5. Firstly, in a process 21, a steady 
functional model is obtained which reproduces a function, performance, 
and characteristic of a steady state from the functional model of the 
25 part. This steady functional model turns into a model including only a 
steady characteristic except a transient characteristic. In a process 22, 
a steady test of a part which forms an identification subject is 
performed to collect steady test data. In a process 23, the steady 
internal characteristic value is identified by using the test data to be 
30 provided to the functional model by a process 24. 

The transient identification process 3 includes the following 
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processes according to claims 6-9. Firstly in a process 31, a transient 
test model of a part which forms an identification subject is prepared, 
and the test is performed to collect transient test data. In a process 32, 
a transient internal characteristic value of the functional model is 
5 identified by using the test data. As shown in a process 4, the steady 
internal characteristic value identified in the process 24 is reflected to 
the transient identification process 3. 

Hereinafter, along these processes, the test model and the 
identification method for identifying the functional model will be 
10 described taking a motor as an example enabling a brief description, 
ju Firstly, the functional model of the process 1 common to each 

2 claim will be described by taking a DC (direct current) motor as an 

example. 

£0 

SI <Modeling of motor> 

m 

15 (1) Basic function of motor 
L A basic function of a motor is that electric energy is inputted, and 

W the energy is converted into a rotation energy to be outputted. As for 

M= the energy conversion at this time, a voltage V M and a current I M are 

U applied to the electric system of the motor whereby an angular velocity 

20 co M and a torque T M are outputted from the rotation system. Also, a 
loss characteristic, an accumulation characteristic, and an additional 
load are included in the motor as an internal load, so that these 
determine a behavior and a loss characteristic of the motor. 

A torque coefficient x T [Nm/A] expressed by the ratio of current 
25 and torque, and a velocity coefficient % „ [V/(rad/sec)] expressed by the 
ratio of angular velocity and induced voltage are known between the 
electric system and the rotation system of the DC motor. An ideal 
motor model which directly converts an input/output state value with 
these coefficients being made a basic function can be expressed by the 
30 following equation: 
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In equation (1), P is the number of poles in armature, "a" is the 
number of armature parallel circuits, Z is the number of all conductors, 
and 0 is all of the magnetic flux per pole. They are characteristics 
5 determined by a motor structure. Also, the magnetic flux <2> is 
determined by the structure in a separately excited type motor which 
uses a permanent magnet for a field in the same way as above, and is 
influenced by a state value of the electric system in a self-excited motor 
having a field winding. Eq.(l) is a basic one common to a DC rotator 
10 applied to a DC generator as well. 

A torque coefficient % T and a velocity coefficient % a of Eq.(l) 
4] performs a physical unit system transformation from current into 

5 torque, or from angular velocity into voltage in the same equation, 

fy Therefore, these coefficients will be hereafter called motor constants 

LI 15 expressed by M M . As for the relationship between the motor constants 
P and the internal load, the electric system and the rotation system are 

considered to be the motor internal load with the coefficients being 
bordered. A basic form of the motor model combining the relationship 
of the basic function of the motor constant M M and the internal load is 
20 shown in Fig. 2. 

In Fig.2, the input/output state value of the electric system is 
connected to a power supply, and the input/output state value of the 
mechanical system is united with an external machine load. Also, as 
for the flow of the state value, the lower flow represents a flow system 
25 in which the current I M of the power supply is converted by M M into a 
flow of the output torque T M which drives the external load, while the 
upper flow represents a potential system in which the angular velocity 
co M received from the external load is converted into the motor induced 



5 



voltage by M M to turn into a flow of the voltage V M returned to the 
power supply. 

Thus, the mechanical system and the electric system are united 
by the strength of energy called a potential quantity and the quantity 
5 of energy called a flow quantity, thereby enabling energy rules of two 
dimensions, equilibrium of forces and the continuity of velocity, to be 
satisfied. 

(2) Functional model of motor 

An electric circuit diagram of the motor model is shown in Fig.3. 
10 In the electric system of Fig.3, R M is a winding resistance, Lm is 

tt an inductance, and R c is an insulation resistance of the winding. Also, 

5 the rotation system has a moment of inertia J M , a viscous resistance 

£ D M , and a friction torque Tmf. The functional model expressing these 

3 relationships is shown in Figs.4(l) and 4(2). It is to be noted that such 

^15 a functional model itself has already been disclosed in the Japanese 
jf. Patent Application Laid-open No. 9-91334. 

In Fig.4(l), the right side of the motor constant M M forms the 
rotation system internal load, and the left side thereof forms the 
electric system internal load. This model has a voltage drop E MB of a 
20 brush which supplies the current I M to the armature, and a friction 
torque T MF of the rotation system as the additional load. When the 
influence of both is neglected, Tmf and Emb can be assumed to be 0. It is 
to be noted that x M and x L in Fig. 4(1) are internal state values. Also, 
the following velocity correction coefficient 6 M is added to the motor 
25 constant M M by taking into account the variation of the torque 
coefficient % T and the velocity coefficient x 
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A government equation of such a motor functional model in 
Fig. 4(1) is given by the following equation: 
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In equation (3), the lst-2nd rows form state equations, and the 
2nd-4th rows form output equations. 

It is to be noted that Eq.(3) can be derived from the following 
5 process shown in Fig.4(l); 

- A white circle adds the inputted state values to make the output 
state value; 

S3 

5 • At a black circle, the outputted state value is equal to the inputted 

m state value; 

7i 10 • A characteristic expressed by a square makes the product of the 

* input state value and the characteristic the output state value; 

M= • A white circle with a mark X outputs the product of the input state 

ru 

m values or input signals; 

K • A while circle with a mark - makes the state value negative; 

15 - A triangle integrates a differential quantity of the inputted state 
value to make the output state value; 

• A large white circle on which a variable name is written expresses 
the state value generated in the functional model as an additional load. 
(3) Uniting with motor power supply 

20 A battery forming the power supply of the motor is modeled. Fig. 5 

shows a circuit diagram of the battery, and Fig.6 shows its functional 
model. In the battery of Fig.6, a voltage V B and a current I B are 
connected to the voltage V M and the current I M of the motor shown in 
Fig.4(l). Also, the characteristic of Fig.6 has E 0 as an induced 

25 electromotive force of the battery, and R E as an internal resistance of 
the battery. 

The government equation of the battery in Fig. 5 is given as 



7 



follows: 
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Uniting the motor model and the battery model is given by the 
following equation by substituting equation (4) for Eq.(3). 



•Eq.(5) 



The lst-2nd rows of equation (5) form the state equations, and 
the 3rd-5th rows form the output equations. An abstracted equivalent 
admittance Y 0 and a current source I 0 in Eq.(5) are expressed by the 
following equations: 
1 
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Eq.(6) 



^0 = Y Q E 0 



•Eq.(7) 

Thus, a characteristic value identification method according to 
the present invention uses a functional model of a part based on a 
potential quantity and a flow quantity representing energy applied to 
the part, comprises a first process for preparing that functional model, 
a second process for converting the functional model into a steady 
functional model in a steady state to identify a steady internal 
characteristic value, and a third process for identifying a transient 
internal characteristic value in a transient state of the functional 
model by using the steady internal characteristic value (claim 1). 

Thus, in the characteristic value identification method of the 
present invention, the functional model of the part is prepared by the 
energy rules called the potential quantity and the flow quantity, and 
the transient identification is performed after the steady identification 
of its internal characteristic value, thereby enabling accurate 
identification of the characteristic values for all the parts and modeling 
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of an integrated product. 

The above-mentioned second process may include a first step for 
determining an internal characteristic value of at least one steady test 
model from the steady functional model, a second step for collecting 
5 steady test data by performing a test corresponding to the steady test 
model, and a third step for identifying a steady internal characteristic 
value of the internal characteristic value based on the steady test data 
(claim 2). 

The above-mentioned first step may determine the internal 
10 characteristic value from a government equation in the steady state of 
the functional model (claim 3). 

Also, the above-mentioned third step may convert the 
G) government equation into a recurrence equation to determine the 

LH steady internal characteristic value from a recurrence coefficient of the 

T 15 recurrence equation (claim 4), and may divide the steady internal 
Sf s characteristic value into a known factor and an unknown factor to 

W identify the steady internal characteristic value of the unknown factor 

O (claim 5). 

Furthermore, the above-mentioned third process may include a 
20 first step for determining an internal characteristic value of at least 
one transient test model in a transient state of the functional model, a 
second step for collecting transient test data by performing a test 
corresponding to the transient test model, a third step for applying the 
steady internal characteristic value to the internal characteristic value 
25 of the transient test model to generate transient phenomenon 
reproduction data, and a fourth step for correcting the transient 
phenomenon reproduction data based on an error between the 
transient phenomenon reproduction data and the transient test data, 
thereby identifying a transient internal characteristic value (claim 6). 
30 When the error does not he within an allowable range, the 

above-mentioned fourth step may repeatedly correct a predetermined 
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transient internal characteristic value within the transient 
phenomenon reproduction data until the error lies within the 
allowable range, and determine the transient internal characteristic 
value to be identified when the error lies within the allowable range 
5 (claim 7). 

Also, the above-mentioned fourth step may preliminarily 
calculate a variance deviation, as a time history sensitivity, to an 
initial value at a time when each transient internal characteristic 
value is increased or decreased at a fixed ratio, and select a transient 
10 internal characteristic value having a maximum sensitivity within the 
^ time history sensitivity as the predetermined transient internal 

° characteristic value (claim 8), or may select a transient internal 

O 

y* characteristic value having the time history sensitivity similar to the 

sj error as the predetermined transient internal characteristic value 

]p 15 (claim 9). 

!\ It is to be noted that the fourth step may simultaneously select a 

PJ plurality of transient internal characteristic values having different 

P maximum sensitivity times as the predetermined transient internal 

characteristic value (claim 10). 
20 Also, a characteristic value identification apparatus according to 

the present invention comprises block replacement means for a 
functional model of a part, test reproduction means for reproducing at 
least one steady test model in a steady state of the functional model 
and at least one transient test model in a transient state, testing 
25 means of the part for performing a steady test and a transient test 
respectively corresponding to the steady test model and the transient 
test model, measurement means for collecting steady test data and 
transient test data at a time when a steady test and a transient test of 
the part are performed by the testing means, and calculating means for 
30 identifying a steady internal characteristic value of the steady test 
model by using the steady test data, for applying the steady internal 
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characteristic value to the transient test model to generate transient 
phenomenon reproduction data, and for correcting the transient 
phenomenon reproduction data based on an error between the 
transient phenomenon reproduction data and the transient test data, 
thereby identifying a transient internal characteristic value (claim 11). 

Thus, it becomes possible to quickly identify the functional model 
of the same kind. 

When the error does not lie within an allowable range, the 
above-mentioned calculating means may repeatedly correct a 
predetermined transient internal characteristic value within the 
transient phenomenon reproduction data until the error lies within the 
allowable range, and determine the transient internal characteristic 
value to be identified when the error lies within the allowable range 
(claim 12). 

Also, the above-mentioned calculating means may preliminarily 
calculate a variance deviation, as a time history sensitivity, to an 
initial value at a time when each transient internal characteristic 
value is increased or decreased at a fixed ratio, and select a transient 
internal characteristic value having a maximum sensitivity within the 
time history sensitivity as the predetermined transient internal 
characteristic value (claim 13), or may select a transient internal 
characteristic value having the time history sensitivity similar to the 
error as the predetermined transient internal characteristic value 
(claim 14). 

It is to be noted that the calculating means may simultaneously 
select a plurality of transient internal characteristic values having a 
different maximum sensitivity time as the predetermined transient 
internal characteristic value (claim 15). 

Furthermore, a virtual testing system according to the present 
invention incorporates the functional model, as a virtual prototype, 
having an internal characteristic value identified by a characteristic 

11 



value identification apparatus claimed in claim 11 and comprises 
condition assigning means for assigning a driving operation condition 
and an environment condition to the characteristic value identification 
apparatus, observation means for observing reproduction data 
5 obtained by the virtual prototype when the driving operation condition 
and the environment condition are assigned, and evaluation means for 
evaluating an observation result of the observation means (claim 16). 

Namely, in the identification of the internal characteristic value 
by the present invention, an evaluation of a function, a performance, 
10 and a characteristic of a part as well as a product in an actual machine 

El test can be regarded as a virtual test performed on a computer, and the 

O 

D actual machine test performed in a development process of designing, 

b prototyping, and testing in the prior art can be omitted, so that it 

yi becomes possible to shorten the period and to reduce the development 

* 15 cost by the virtual test. 

H The above-mentioned virtual testing system may further 

fU comprise another measurement means for measuring actual machine 

S test data at a time when the driving operation condition and the 

environment condition are provided to an actual machine which forms 
20 a subject of the virtual prototype, and re-identification means of the 
virtual prototype, wherein the evaluation means may compare an 
output of the measurement unit and the observation result, and make 
the re-identification means re-identify the virtual prototype according 
to the comparison result (claim 17). 
25 A fixed virtual prototype may be incorporated into a part of a 

drive system and a load system connected to the part as the above- 
mentioned virtual prototype, the testing means may perform a test 
corresponding to the fixed virtual prototype, and the evaluation means 
at this time may make the re-identification means perform a re- 
30 identification according to the comparison result (claim 18). 
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Brief Description of Drawings 

Fig.l is a block diagram showing a basic concept of a 
characteristic value identification method and an apparatus therefor 
according to the present invention; 
5 Fig. 2 is a block diagrams showing a basic system of a motor 

model with a potential quantity and a flow quantity of energy being 
mutually related; 

Fig. 3 is a diagram showing a general electric circuit of a motor; 
Fig. 4 is a diagram respectively showing a functional model and 
10 symbols of a motor used in a characteristic value identification method 
and an apparatus therefor according to the present invention; 

Fig. 5 is a diagram showing a general electric circuit of a battery; 
"Tig. 6 is a diagram showing a functional model of a battery; 
m Fig.7 is a block diagram showing a process for preparing a steady 

y 15 test model group from a motor steady functional model relating to a 
characteristic value identification method and an apparatus therefor 
according to the present invention; 

Fig. 8 is a diagram showing a steady functional model of a motor; 
Fig. 9 is a diagram showing a motor model at the time of a steady 
20 load test; 

Fig. 10 is a diagram showing a motor model at the time of a load 

test; 

Fig. 11 is a diagram showing a motor model at the time of a lock 

test; 

25 Fig. 12 is a diagram showing a motor model at the time of a 

induced voltage test; 

Fig. 13 is a diagram showing a motor model at the time of a 
braking test; 

Fig. 14 is a block diagram showing an identification process of a 
30 steady functional model from a steady test model group in a 
characteristic value identification method and an apparatus therefor 
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according to the present invention; 

Fig. 15 is a graph showing a motor characteristic of a steady state; 
Fig. 16 is a diagram showing an external shape of a motor; 
Fig. 17 is a graph showing a steady internal characteristic value 
5 of a motor with an ordinate and an abscissa in Fig. 15 being exchanged; 

Fig. 18 is a graph showing a steady internal characteristic value 
of a motor; 

Fig. 19 is a block diagram showing a transient identification 
process of a functional model in a characteristic value identification 
10 method and an apparatus therefor according to the present invention; 

N> Fig. 20 is a diagram showing a structure of an electrically driven 

□ 

0 arm mechanism; 

S Fig. 21 is a diagram showing a model at the time of an actual load 

^ transient test of a motor; 

yi 

1 15 Fig.22 is a diagram showing a transient test model at the time 
when a motor axis is fixed; 

Fig.23 is a diagram showing a fixed torque transient test model of 
a motor; 

Fig.24 is a diagram showing a model at the time of an inertial 
20 rotation test of a motor; 

Fig.25 is a diagram showing a motor operation switch circuit; 
Fig. 26 is a diagram showing a functional model of a motor 
operation switch; 

Fig.27 is a diagram showing an electric circuit of a switching 

25 method; 

Fig.28 is a diagram showing a functional model of an electric 
circuit of a switching method; 

Fig.29 is a diagram showing an electric circuit of a fixed current 

control; 

30 Fig.30 is a diagram showing a functional model of a fixed current 

control circuit; 
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Fig. 31 is a graph showing a trial simulation result before 
performing a transient identification in a characteristic value 
identification method and an apparatus therefor according to the 
present invention; 

5 Fig. 32 is a flow chart illustrating an operation of a transient 

identification in a characteristic value identification method and an 
apparatus therefor according to the present invention; 

Fig. 33 is a graph respectively showing a sensitivity of a 
characteristic value in a characteristic value identification method and 
10 an apparatus therefor according to the present invention; 

Fig. 34 is a graph respectively showing a relationship between a 
current deviation and a characteristic value time history sensitivity in 
a characteristic value identification method and an apparatus therefor 
according to the present invention; 
15 Fig.35 is a graph respectively showing a simulation result after 

performing an identification in a characteristic value identification 
method and an apparatus therefor according to the present invention; 

Fig. 36 is a graph respectively showing a simulation result of a 
start and a stop in a characteristic value identification method and an 
20 apparatus therefor according to the present invention; 

Fig. 37 is a graph respectively showing a simulation result of a 
regular/reversed rotation and brake in a characteristic value 
identification method and an apparatus therefor according to the 
present invention; 

25 Fig. 38 is a block diagram illustrating a characteristic value 

identification apparatus according to the present invention by taking a 

motor as an example; 

Fig. 39 is a block diagram showing a relationship between a 

product development and a model; 
30 Fig. 40 is a block diagram showing a concept of a virtual test 

according to the present invention; 
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Fig.41 is a diagram showing a model at the time when a virtual 
testing system according to the present invention is applied to a car; 

Fig. 42 is a block diagram showing an embodiment of a virtual 
testing system according to the present invention; 
5 Fig.43 is a block diagram showing an example applying a virtual 

testing system according to the present invention to an actual machine 
test of a transmission; 

Fig. 44 is a model diagram showing a virtual testing system of a 
drive system, a transmission, and a load system; and 
10 Fig.45 is a block diagram showing an example applying a virtual 

testing system according to the present invention to an actual machine 
test of a power train. 
{J Description of Reference Numerals 

M 1 motor model (block replacement structure) 

5 15 2 test reproduction model 
L 3 tested motor 

4 drive circuit unit 
= 5 load generator 

6 electric load 
20 8, 34 measurement unit 

9 calculation unit 

10 virtual test model 

11 virtual prototype 

12 test standard model 

25 13 engine model (motor model) 
14 body model (generator model) 

21 drive model (battery operation model) 

22 load model (load generator model) 

23, 37, 40 observation model, observation system 
30 24 characteristic value update model 
25 model replacement unit 
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30 virtual test portion 

31, 38, 39 drive operation model, operation system 

32, 41 environment condition model, environment system 
33 actual machine test portion 

5 35 evaluation model 

36 control system input/output unit 

50 tested transmission 

51 tested power train (P/T) 
71 drive system control unit 

10 72 load system control unit 

Throughout the figures, like reference numerals indicate like or 
corresponding components. 



rj Best Mode for Carrying Out the Invention 

42 15 Hereinafter, an embodiment of a characteristic value (model) 

s identification method and an apparatus according to the present 

invention schematically shown in Fig.l will be described. 

rf l Functio nal model in steady state 

p In order to perform the steady identification process 2 shown in 

20 Fig.l, the steady functional model conversion is required as shown in 
the process 21. In this steady functional model conversion, it is 
necessary as will be described hereafter to convert the functional 
model into the steady functional model for reproducing the steady 
state of the part, and to further convert the steady functional model 
25 into the steady test model according to steady tests of the part (process 
210). 

Since the steady identification utilizes the steady test data, the 
steady test model can make the steady functional model of the 
identification subject accord with the tested object to which the actual 
30 machine test is performed, by modeling the tested object according to 
the condition and the contents of the steady test for generating the test 
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data. Fig.7 shows a flow of a series of model conversions. 

1. 1 Conversion from f unction model into steady functional model 

Since there is no transient variation in a steady state of a motor, 
neither the inductance Lm nor the moment of inertia J M shown in Fig. 
4A influences the steady internal characteristic value. Also, since the 
insulation resistance R c in Fig. 4A has an extremely large value, the 
influence to the motor current I M can be neglected. From the above- 
mentioned points, the motor model in the steady state where Lm, J m , R c 
of Fig.4(l) are eliminated turns into a functional model shown in Fig.8. 

The left side of Fig.8 shows a battery model of Fig. 6 in which the 
operation switch of the electric system is omitted, and the right side 
thereof shows a motor model in a steady state. From Fig.8, the motor 
model in the steady state can be expressed by the following 
government equation: 
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•Eq.(8) 



For another method of deriving the steady government equation, 
the government equation in the steady state obtained by eliminating 
Lm, Jm, Rc &° m the government Eq.(3) of the functional model in 
Fig.4(l) can be obtained by the following processes: 

® The equation obtained by eliminating the differential state values 
of x M and x L from Eq.(3) with Lm and J M being assumed to be 0 is 
given as follows: 
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(3) The equation obtained by eliminating the internal state values of 
25 x M and x L from equation (9) is given as follows: 
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• -Eq.(lO) 

"D" in equation (10) is assumed to be the following equation: 
D = D M (R M + /O + ^M w 2 • * • Eq. (11) 

5 ® The equation obtained by eliminating the insulation resistance by 
making R c =00 from Eq.(10) is given as follows: 
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The government equation (12) in the steady state derived from 
the government Eq.(3) of the functional model is equal to the 
government equation (8) derived from the functional model in the 
steady state in Fig. 8. 

Accordingly, there are two methods for obtaining the steady 
government equation reproducing the steady state of the part; one by 
converting the functional model reproducing the transient state of the 
part into the steady functional model, and the other by deriving the 
steady government equation from the government equation of the 
functional model. This indicates that the functional model in the 
steady state can be reversely depicted from the derived steady 
government equation, and that there is a reversibility between the 
steady functional model and the steady government equation. 
1 2 Oonypvftirm from steady f u nctional model into test mode}, 

A steady state of an actuality by the steady functional model will 
be considered, and the test model will be examined. The test model in 
the steady state models the test of the actuality which provides a 
specified condition to an input/output system of a tested object, from 
the steady functional model and its government equation. The steady 
internal characteristic value is identified by relating the test data 



19 



obtained when the specified condition is assigned to the tested object, 
with the government equation. 

The relationship between the specified condition in which the 
steady internal characteristic value is assigned to the motor model and 
5 the test is as follows: 

© In case a fixed torque is provided to the motor with the power 
supply being connected, the model reproduces the relationship between 
the load applied to the motor and the steady internal characteristic 
value, so that the relationship between the motor current, the angular 
10 acceleration, and the load torque can be observed (steady load 
g test:Fig.9). 

P (D The case where the input torque T M is assumed to be 0 corresponds 

60 to a motor rotation in an idle state, so that an idle angular velocity and 

m an idle current are observed (idle test:Fig. 10). 

f 15 (D The case where an output angular velocity co M is assumed to be 0 
!f« corresponds to a lock state or a state at a starting moment when the 

W motor rotation is compulsorily stopped, so that a lock current and a 

=5= _____ 

b lock torque can be observed (motor lock test:Fig.ll). 

- . 

© In case the input current I M is assumed to be 0 and the input torque 
20 T M is provided, the model which reproduces a generator function of the 
motor is obtained, so that the induced voltage and the angular 
acceleration can be observed (induced voltage test:Fig.l2). 
(D In case a resistor is connected to the input/output state value of the 
electric system to provide a fixed torque, the model which reproduces a 
25 regenerative braking force of the motor can be obtained, so that the 
relationship between the motor current, the angular acceleration, and 
a braking torque can be observed (steady braking test:Fig.l3). 

These test models are different from each other depending on 
how to provide the specified conditions. However, in an actual 
30 modeling, the values of the steady characteristic such as a steady 
functional model or a steady load test model and the input state are 
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made 0 or an infinite value (large value within a range where the value 
does not influence a reproduction result), so that each test model can 
be equivalently realized. It is necessary to perform at least one of such 
steady tests. 

In view of the above-mentioned description, the motor test model 
performed from the steady functional model in Fig.8 under the 
specified conditions, is given as follows: 
1 2 1 Motor steady load t est 

Fig.9 shows a steady load test which observes the performance 
and the characteristic by providing the load torque to the motor. 

In Fig.9, the input torque T M is provided as a drive torque by the 
additional load from the steady functional model in Fig.8, and the 
output angular velocity o> M is represented as an observation quantity. 
Accordingly, the government Eq.(8) of the motor is transformed into 
the following equation: 



The government equation of the steady load test obtained by 
mutually substituting the government Eq.(4) of the battery for 
equation (13) is given as follows: 




••Eq.(13) 



1 MuiEt-E^-iRs+RHWH+TMr) 




•Eq.(14) 
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Equation (14) is a government equation reproducing a rated load 
test of the motor, where all of the equations assume the output state 
value reproducing the state value of the steady state. 
1 . 2, 2 M otor idl e test 

In the motor idle state, the input torque T M of the steady load test 
model in Fig.9 becomes 0. The idle test model where the input torque 
T M is eliminated from Fig.9 can be shown by Fig. 10. 

The government equation of the motor in Fig. 10 turns into the 
following equation obtained by eliminating the load torque T M from the 
government Eq.(13) of the steady load test: 



L M 



Eq.(15) 



The government equation of the idle test obtained by mutually 
substituting the government Eq.(4) of the battery for equation (15) is 
given as follows: 



1 M M < E 0 -E^ ) - {R B +R M 



£ 0 + 1 A 



R. 



B 



A, 



MF 



1 + 1 ^ + D M R B 



— Imf 



A, 



Eq.(16) 



L 2. 3 Motor lock test 

At the starting moment of the motor, the angular velocity co M 
becomes 0, where it can be regarded as a motor lock state in which the 
motor rotation is compulsorily stopped. Accordingly, the functional 
model of the lock test reproducing the motor lock can be shown by 
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Fig. 11 where the output angular velocity o> M of the steady functional 
model in Fig. 8 is supposed to be 0. It is to be noted that since a brush 
voltage drop E MB generated in the stop state of the motor rotation does 
not generate a fixed voltage drop but only assumes a minute voltage 
drop by a brush resistance, it is neglected and omitted. 

In Fig. 11, each steady internal characteristic value relating to the 
output angular velocity a> M is eliminated from the steady functional 
model in Fig. 8, so that the input torque T M is shown by the observation 
quantity of the lock torque. The government equation of Fig. 11 is given 
as follows: 

" [/J -Eq.(17) 

The government equation of the lock test obtained by 
substituting the government Eq.(4) of the battery for equation (17) to 
be integrated is given as follows: 









J*. 







-M 



M 



[1] 



■Eq.(18) 



1 2. 4 Moto r induced voltage test 

Fig. 12 shows a model of an induced voltage test which provides 
the torque to the motor to observe a counter electromotive force V M 
induced by the armature and the angular velocity co M 

In Fig. 12, the steady internal characteristic value relating to the 
input current I M is eliminated from the steady functional model in 
Fig.8, so that the output angular velocity co M is shown by the 
observation quantity. The government equation of Fig. 12 is given as 
follows: 
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•Eq.(19) 



1 ?. K Motor steady braking test 

Fig. 13 shows a steady model of a regenerative braking which 
provides the torque to the motor and reversely returns the counter 
electromotive force V M induced by the armature to the motor through a 
braking resistance R L to provide a braking force. 

In Fig. 13, the battery is eliminated from the steady functional 
model in Fig.8, the braking resistance R L is connected, and the input 
torque T M of the motor is provided as a drive torque. In addition, the 
output angular velocity a) M is shown as the observation quantity. The 
government equation of the motor in Fig. 13 is given as follows: 



\Jm J 



(R L +R M ) + 



M 



A. 



Rr 



D. 







[1] 



*Eq.(20) 



The above-mentioned test models can be equivalently modeled by 
setting the following values to the characteristic or the input state 
value of the steady load test model shown in Fig.9: 

® The idle test model becomes equivalent by rendering the input 
torque T M =0. 

(D The lock test model becomes equivalent by rendering a viscous 
damping coefficient and the brush voltage drop respectively D M =<*>, 

(B) The induced voltage test model becomes equivalent by rendering 
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the motor current and the brush voltage drop respectively I M =0, E MB =0. 
(D The braking test model becomes equivalent by rendering the 
internal electromotive force of the battery E 0 =0, and by rendering the 
internal resistance R B the braking resistance R L . 
5 2- Identification of steady sState 

Fig. 14 shows an identification process 2 of the steady internal 
characteristic value included by the steady functional model having 
been described referring to Fig.l. Since the steady internal 
characteristic value known beforehand from a design value, the 
10 characteristic value of a similar article, the measured value of the 
u characteristic, and the like can be determined as a fixed value in the 

5 steady identification process (process 26), remaining unfixed values 

t can be identified (process 27). 

SI As for an analysis technique applied to this, the government 

£ 15 equation of the steady test model is converted into a statistical model 
L such as a canonical correlation analysis or a multiple fold recurrence 

H analysis (simple recurrence analysis), determines a deviated 

M= recurrence coefficient from the steady test data, and derives the steady 

p internal characteristic value from the relationship between the 

20 deviated recurrence coefficient and the government equation (process 
28). Also, as for the steady characteristic formed by many factors, the 
factors are analyzed, and divided into factors of which quantification is 
possible and factors of which quantification is difficult by modeling the 
forming process of the steady internal characteristic value with a 
25 theoretical equation, an experimental equation, and the like, thereby 
identifying the factors of which quantification is difficult (process 29). 
9. 1 TRst plan of steady identification 

In the steady functional model shown in Fig. 8, there are six kinds 
of the steady internal characteristic values M M , R M , D M , 5 m, TW, and 
30 Emb which form the subject of the steady identification. In order to 
determine these steady internal characteristic values, it is necessary to 
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firstly clarify the function of the steady model in Fig.8 and the 
characteristic of the steady internal characteristic values, and to 
clarify the relationship with the test method. 

9 t 1 , 1 VftrificRtion of kn o wn/mpasurahle steady internal characteristic 
5 y?hift (process 26) 

Firstly, known characteristic values and the characteristic values 
which can be solely measured are clarified in the steady internal 
characteristic values included in the steady functional model. Among 
the characteristic values which govern the steady state of a general 
10 product or part, there are a large number of characteristic values 
u which can be clarified without performing the steady test such as a 

2 characteristic value clarified in the literature or the like, a 

characteristic value which can directly utilize the design value, and a 
Si characteristic value which can be directly measured. Accordingly, it is 

% 15 important to determine the characteristic values by a preliminary 
J\ investigation or measurement in case of the steady identification. 

fU Reviewing the motor used here, such characteristic values to be 

S determined are a brush voltage drop Emb, an insulation resistance R c 

S omitted at the time of modeling the steady functional model in Fig.8, a 

20 motor winding resistance R M , and the like. These steady 
characteristics are given as follows: 

© Since being a voltage drop due to a graphite brush, the brush 
voltage drop E MB assumes constant regardless of the current I M during 
the motor rotation, and assumes a minute DC resistance in the stop 
25 state. Accordingly, 0.7 [V] shown in the literature is applied during the 
motor rotation, and the direct current in the stop state can be 
neglected. 

(2) The insulation resistance R c assumes the value [MQ] which does 
not influence the result. 
30 (3) The winding resistance R M of the motor can be determined by the 
ohmmeter or by calculating the motor voltage and current in the motor 
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lock state shown in Fig. 11. 

9 1 2 Steady f unctional mndRl and statistical model (process 28) 

There is a multivariate analysis such as a multiple fold 
recurrence analysis and a simple recurrence analysis in the analysis 
5 technique frequently used for the product development. The recurrence 
equation of the multiple fold recurrence model applied thereto is given 
as follows: 

y = b 0 + Vi + Mi— M« + e rr -Eq.(21) 
Applying the multiple fold recurrence equation shown in equation 
10 (21) to the steady functional model including a plurality of unknown 
steady internal characteristic values will be examined. As for the 
method of application, dependent variables and independent variables 
in the steady government Eq.(8) are made correspond to an objective 
variable "y" and predictor variables x r x n , to be determined from the 
i 15 steady test data. 

Also, deviated recurrence coefficients b 0 -b n are values combining 
known and unknown steady internal characteristic values. Accordingly, 
unknown steady internal characteristic values can be derived by 
relating (associating) the deviated recurrence coefficients to (with) the 
20 steady internal characteristic values. In Eq.(21), e rr is a residual. Also, 
when the multiple fold recurrence equation has a plurality of different 
objective variables, these are expressed by simultaneous equations so 
that the canonical correlation analysis is applied. As a matter of course, 
the test data obtained by a plurality of different test methods, the 
25 deviated recurrence coefficient, and the steady internal characteristic 
value of the recurrence model can be related to each other. 

In the multiple fold recurrence analysis of the motor, the steady 
internal characteristic value can be clarified by the following process 
from the data of the steady load test shown in Fig.9: 
30 Firstly, as for the multiple fold recurrence model, the following 

multiple fold recurrence equations can be derived from the government 
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Eq.(13) of the steady load test model in Fig.9: 



— Eq.(22) 




-Eq.(23) 



It is to be noted that equation (22) corresponds to the 1st row of 
Eq.(13), and equation (23) corresponds to the 2nd row of Eq.(13). Since 
the internal resistance R B of the battery is very small, and the internal 
electromotive force E 0 of the battery and a motor output voltage V M are 
considered to be approximately equal, E 0 in Eq.(23) is substituted for 
the dependent variable V M (motor output voltage) of Eq.(13). 

The recurrence equation with the torque T M of these equations 
being made an objective variable and the motor current I M and the 
angular velocity co M being made predictor variables is given as follows 
(process 281), where the upper equation of the following equations (24) 
forms a multiple fold recurrence equation with two predictor variables, 
and the lower equation forms a simple recurrence equation with a 
single predictor variable: 



Also, constants (intercepts) B 10 , B 20 , and deviated recurrence 
coefficients B u , B 12 , B 21 can be expressed by the following equations: 



The steady internal characteristic values derived from the 
deviated recurrence coefficients shown in equations (25) are given as 
follows (process 282): 

Firstly, Tmf, D m , M m correspond to the respective deviated 
recurrence coefficients from the upper four rows of Eqs. (25). 




•Eq.(24) 




MF 



•••Eq.(25) 



Ml 



T MF — B ]0 
D M =-B U 
M M =B n 



Eq.(26) 



Then, 6 M can be derived by substituting equation (26) for four 
rows of Eqs. (25) as follows: 

5 »= TB B B W (£ °- g -> -Eq.(27) 

It is to be noted that the induced voltage E 0 of the battery and the 
brush voltage drop E MB in equation (27) are assumed to be known 
steady internal characteristic values. 

Moreover, R M can be derived from the following equation by 
substituting Eqs. (26) and (27) for five rows of Eqs. (25): 



S 10 R M =%^(E 0 -E^) -Eq.(28) 



While in the above-mentioned example, all of the steady internal 
characteristic values have been determined from the government 
Eq.(8) by the recurrence analysis, in many cases of the actual steady 
identification, the steady internal characteristic values are derived by 

15 combining a plurality of test methods. For example, in the above- 
mentioned example, as for the 3rd row of Eqs.(26) deriving the motor 
constant M M and equation (28) determining the winding resistance R M , 
both M M and R M can be determined from Eq.(17) of the lock test model 
and the test data as well. 

20 2. 1 . 3 Factor analysis of steady characteristic (process 29) 

The steady internal characteristic value of the steady functional 
model includes many internal factors forming this characteristic. 
Accordingly, it is necessary to clarify the factors composing each 
characteristic value by performing a factor analysis, and to formulate 

25 the result. Also, it is necessary to divide the factors into known factors 
and unknown factors based on the result, and to perform a steady 
identification for clarifying an influence degree of the unknown or 
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uncertain factors by the model or the test data. 

To the steady identification method, a sensitivity analysis or the 
like is effective by which the relationship between the variance of the 
government factor and the steady internal characteristic value is taken 
5 as a sensitivity. Particularly, in case the steady internal characteristic 
value has a non-linear characteristic influenced by another state value, 
the sensitivity analysis is important. 

If the motor constant is considered, for example, P, "a", and Z are 
known values in Eq.(l) since they are design constants determined by 
10 the motor structure. However, it is difficult to determine all of the 
magnetic fluxes 0 per pole since it is distributed in the gap between 
the field and the armature. Accordingly, to identify the motor constant 
means to identify the magnetic flux0 by the test data. All of the 
5 magnetic fluxes 4> per pole determined from the identified motor 

W 15 constant M M are given as follows: 

^ =Mjw || x l0 8 -Eq.(29) 

Moreover, as a brief example, the winding resistance R M of the 
motor is given by the following equation determined by a line length T, 
a sectional area "S", a specific resistance " p and a temperature 
20 coefficient " a being governed by a temperature "t": 

R M = P Ul+ca) -Eq.<30) 

From equation (30), 1, S, p , and a assume the design values 
determined by the structure and the stuff, while the temperature "t" is 
governed by the temperature of the test environment and the 
25 temperature rise by self-heating of the motor. Accordingly in the test 
for performing the identification, the temperature of the motor forms 
an important factor which governs the characteristic value. Also, upon 
the examination of the high temperature environment and low 
temperature environment, it is required to perform a temperature 
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correction, by Eq.(30), to R M in the government Eq.(3) and the steady 
government equations of equation (8) and subsequent equations. 
2 1 4 Motor steady characte ristic 

The steady internal characteristic value of the motor is 
5 represented by a torque -current characteristic representing the 
relationship between the output torque T M [Nm] and the motor current 
I M [A], and by a torque-velocity characteristic representing the 
relationship with the angular velocity co M [rad/sec], based on the 
result measured in the above-mentioned steady load test (Fig.9). 
10 Firstly, the torque-velocity characteristic can be expressed by the 
following equation obtained by substituting the motor current I M of 
Eqs.(22) and (23), the latter itself representing the torque-current 
characteristic: 

M^- £ » w -H o " + ^K • Eq<31) 

tf 5 15 If the torque-current characteristic and the torque-velocity 

Kl characteristic of the motor are available from a catalog, technical 

5 materials, or the like, the steady internal characteristic value may be 

determined from the characteristic diagram. This motor characteristic 

is shown in Fig. 15. 

20 The relationship between the two recurrence lines in Fig. 15 is 

expressed by Eq.(22). In Fig. 15, the intersection of X axis and the 
torque-angular velocity characteristic shows the motor lock state, and 
the intersection of Y axis and the torque-angular velocity characteristic 
shows the idle state. Accordingly, as for the motor constant M M and the 

25 winding resistance R M their approximate values can be determined by 
substituting the motor current I M , the torque T M , and the voltage V M in 
the lock state for Eq.(17). Also, by substituting the result for Eq.(22), 
an approximate value of the viscous resistance coefficient 8 M at the 
time of idling can be determined, whereas the friction torque Tmf is 

30 neglected. 
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9 9. Motor stead y test and steady identifi cation result 

The characteristic value in the steady state is identified from the 
motor functional model shown in Fig.4(l). It is to be noted that the 
motor used in this identification is a small DC motor of a separately 
5 excited type widely used for an electric actuator of a car, and a 
permanent magnet is used for the field. 

The steady test for the steady identification has been performed 
to the idle test, the lock test, and the steady load test. The idle test and 
the lock test have been performed as a test under the specified 
10 conditions of the motor. The former has been performed with the 
output axis of the motor being released in the state of the input torque 
S T M =0, and the latter has been performed with the output axis being 

rf fixed in the state of the output angular velocity a> M =0. 

W Moreover, in the steady load test, a string hung with a weight has 

W 15 been winded onto a drum attached to the motor output axis, whereby 
T the load torque has been calculated from the mass of the weight and 

% the radius of the drum. Also, the test conditions are shown in the 

!jf following Table 3 together with the result of the steady identification. 

Q Fig. 17 shows a result of the steady test performed to the motor of 

20 Fig. 16 and the recurrence equations, where X and Y axes of Fig. 15 are 
exchanged. 

Fig. 17 is a diagram showing the measured values of the torque- 
current characteristic representing the relationship between the 
output torque T M [Nm] in the steady state and the motor current I M [A], 

25 and representing the relationship with the angular velocity co M 
[rad/sec], and its recurrence equation. The characteristic value on the 
X axis represents the motor idle state (T M =0) as mentioned above, and 
the intersection between the torque-velocity characteristic and the Y 
axis shows the lock state or the starting moment of the motor with the 

30 angular velocity <o M =0. The points shown by variable names in Fig. 17 
show a basic characteristic and a rating characteristic (dotted line) in 
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Table 1. Also, the black circle in Fig. 17 shows the measured values. 

Moreover, the basic performance at the time of starting and idling 
obtained in Fig. 17 is shown in Table 1, and a catalog specification 
published by a motor maker is shown in Table 2 as reference data. 

TABLE 1 



C3 

m 



CHARACTERISTIC NAME 


SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


MOTOR IDLE CURRENT 




[A] 


0. 1 


IDLE ANGULAR VELOCITY 




[rad/sec] 


1410. 0 


STARTING TORQUE 




[Nm] 


28.0X10" 3 


STARTING CURRENT 




[A] 


4.0 


TEST POWER ^RNAL ELLC I ROHO 1 1 Vt 

SUPPLY POWER SUPPLY INTERNAL 
RESISTANCE 




[V] 


12,0 


* 8 


[A] 


1. 0X10~ 3 



TABLE 2 



CHARACTERISTIC NAME 


SYMBOL 


UNIT CHARACTERISTIC VALUE CONDITION 


RATED CURRENT 




[A] 


0. 57 MAX 0. 9 


RATED ANGULAR VELOCITY 




[rad/scc] 


1210.0 TOLERANCE ±160.0 


RATED LOAD 




[Nm] 


4.0X10" 3 


RATED VOLTAGE 




[v] 


12.0 



Table 3 shows a result of the steady identification by the test data 
of Fig. 17. 
TABLE 3 



CHARACTERISTIC NAME 



SYMBOL UNIT 



MOTOR STATE 
CHARACTERISTIC 



VISCOUS RESISTANCE n [Nm/(rad/sec)] 

COEFFICIENT M 
KINETIC FRICTION TORQUE T m? lN m J 



CHAR ACTERISTIC VALUE 
0. 65x10-* 
0. 1 x 10" 6 



WINDING RESISTANCE R M [ Q ] 

BRUSH VOLTAGE E BR [V] 

INSULATION RESISTANCE R c [MQ] 



MOTOR CONSTANT M M 
VELOCITY COEFFICIENT CORRECTION 5 M 



TEST CONDITION 
(KNOWN) 



TEST VOLTAGE 
INTERNAL RESISTANCE 
AMBIENT TEMPERATURE 



Z 0 



[Nm/A] 



[V] 

[Q] 

[deg] 



3.0 
0. 7 
10.0 



7. 25X10' 3 
1.0 



12.0 

50. OX10" 3 
25.0 
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Table 3 is a result assembled according to the process of the 
steady identification examined so far. Among the identified steady 
internal characteristic values, E BR has been assumed to be a prescribed 
value determined from the literature. Also, R c is a measurement result 
5 of the insulation testing set. The other steady internal characteristic 
values M M , R m , D m > 6 M , and T MF are determined from the multiple fold 
recurrence analysis and the simple recurrence analysis by the above- 
mentioned Eqs.(24)-(27). Also, since the internal resistance R B of the 
battery is so small that the resistance value of the current 
10 measurement distributary unit has been applied to the internal 
resistance Rb of the battery, 
p It is to be noted that in order to determine five coefficients B 10 , 

S B 20 , B H , B l2 , and B 2U five simultaneous equations are required. T M> I M) 

F! and co M are vectors where equal to or more than five values 

W 15 determined by the characteristic of Fig. 17 exist. 

Fig. 18 shows a simulation result of the thus identified steady 
hj internal characteristic values by using the steady government 

equation obtained by substituting Eq.(4) of the battery model for Eq.(8) 
P of the motor steady state. 

20 Thick solid lines in Fig. 18 show results of the steady state of the 

angular velocity to M and the motor current I M , while thin dotted lines 
show simulation results performed in Eq.(5) by using the transient 
internal characteristic value before identification described later. It is 
seen from Fig. 18 that the results of the steady state coincides with 
25 each other in the steady region where the transient state becomes 
stable. 

3 Identificati on of transient state (process 32) 

Since the steady internal characteristic values of the motor 
shown in the above-mentioned Table 3 are clarified, the moment of 
30 inertia J M and the inductance Lm as transient internal characteristic 
values governing the transient function of the motor are identified next 
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with the measured value of the motor current I M , by the transient 
identification process 32 shown in Fig. 19. 

As an identification method, approximate values are firstly 
provided to the transient internal characteristic values to perform a 
simulation so that initial values are determined (process 36), then the 
time history sensitivities of the transient internal characteristic values 
for the motor current I M are determined (process 38), and the deviation 
(process 37) between the measured value (process 35) of the motor 
current I M and the simulation value assumes the least for the 
sensitivity indexed. 

Moreover, the viscous resistance coefficient 6 M and the kinetic 
friction torque F w already identified are again identified in the 
transient identification for the confirmation. As for the identification, 
the assembly/disassembly and the replacement by a block replacement 
method can be performed to the model of the operation system and the 
load connected to the motor source supply system and the rotation 
system. 

3 1 Conversion fro m functional model into transient test model 
(process 33) 

The test model in the transient state reproduces an actuality test 
providing a specified condition to the input/output system of the tested 
object in the same way as the steady identification by the functional 
model and its government equation. In the transient identification, the 
test data obtained by providing the specified condition to the tested 
object, and the simulation result by the test model are related, so that 
the characteristic value group included in the government equation is 
identified. Also, the test model of the transient identification models 
the drive or the load to be incorporated by the following method: 
(D Input/output model for driving 

The input/output model for driving is a model for supplying 
driving energy such as a power supply and a torque in order to drive 
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the tested object. This model occasionally has an operation system or a 
control system for operating the supplying method of the driving 
energy. For example, the system connected to the electric system of the 
motor corresponds to such a system including a battery of a power 
supply, an operation switch of a motor, a control unit, and the like. As a 
specific example, there is a braking resistance of the regenerative 
braking mentioned above, or the like. 
(§) Input/output model for load 

The input/output model for the load is a model for accumulating 
or consuming the energy such as a characteristic which forms the load 
of the tested object and an actual load. This model occasionally has the 
operation system or the control system for operating how to provide the 
load energy. In the rotation system or the mechanical system of the 
motor, there is a method for connecting a stiffness (spring stiffness), 
the moment of inertia, the viscous resistance coefficient, and the like. 

It is to be noted that in case the characteristic whose value is 
large is connected, it can be considered that the output axis is fixed. 
Also, as an expediency for providing the fixed load, there is a method 
by which the current, the torque, and the like are provided as an 
additional load, and the voltage and the angular velocity which become 
a pair are made an observation quantity. For operating a machine load, 
the load is switched over by simultaneously using models of a clutch, a 
brake, and the like. 

R 1 1 Actual load tra nsient test model of motor (process 34) 

In the transient test model which reproduces the state of the 
practical use, the battery model in Fig.6 is united with the electric 
system of the motor model in Fig. 4(1), and an actual load model or an 
equivalent load model is further united with the rotation system to be 
converted into a model for an execution closed from the outside. 
(1) Actual load model of electrically driven arm mechanism 

Fig.20 shows a load model of an electrically driven arm 
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mechanism having an operation arm whose operation range is 
restricted by 0 ^ at a rack side of a worm gear, as an example of a 
practical load. 

Fig.21 shows a functional model for an execution in which the 
electrically driven arm mechanism in Fig.20 is incorporated in the 
battery and the motor model. The electrically driven arm mechanism 
in Fig.21 has a stiffness Kc of a stopper, a reduction ratio N G of a 
reduction unit, and a contact radius L c of the stopper. Inside the 
electrically driven arm mechanism, there are a presumed rotation 
angle 9 n of the arm and an internal state x c . 

Also, the condition judgment shown by a hexagon frame in Fig.21 
is a model which judges an arm operation range, and operates a switch 
variable S wc . This functional model has an additional load of torque T L 
as an external load. 

The government equations of the electrically driven arm 
mechanism are given as follows: 



0 -N c . 



No 
0 



x c 
1 



•Eq.(32) 



In this government equations (32), the upper equation indicates 
the government equation of the electrically driven arm mechanism 
shown on the right side of Fig.21, and the lower equation indicates the 
presumed rotation angle 0 L of the arm. Also, the model determining 
the rotation range of the arm is given as follows: 

if 0r * th ™ ( S «c = 1) e!se{S wc = 0) -Eq.(33) 

Equation (33) is for a condition judgment Sw S in which Swc=0 
when the presumed rotation angle £ R of the arm lies within the 
rotation range 0 and otherwise S wc Fl- 

If the above-mentioned Eq.(5) where the battery and motor are 
integrated, and the government Eq.(32) of the electrically driven 
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arm mechanism are integrated, the following equation can be 
obtained: 
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-Eq.(34) 

5 In equation (34), the lst-3rd rows form the state equations, and 

the 4th and the following rows form output equations for observing the 
internal transient state value. It is to be noted that Y 0 and I 0 in the 
equations are an equivalent admittance and a current source 
respectively shown in equations (6) and (7). 
10 (2) Axis fixed model of motor 

The model in which the output axis of the motor is fixed can be 
shown by next Fig. 22. This model can be replaced with the arm 
mechanism on the right side of the motor actual load transient test 
model in Fig.21. 

15 Fig.22 is a model in which the model of the electrically driven 

arm mechanism in the actual load transient test model of Fig.21 is 
replaced with stiffness Kc which fixes the motor output axis. The motor 
torque T MF in Fig. 22 reproduces the lock torque in the state where the 
output axis of the motor is fixed. It is to be noted that the government 

20 equation is omitted in this case. 

3 LU Fixftd torque tran sient test model of motor 

There is a fixed torque test for providing a fixed torque as a load 
in the load test performed in the actual machine test. This test model 
unites the battery model in Fig.6 with the motor model in Fig.4(l) to be 

25 converted into the model for execution closed from the outside. The 
converted functional model for the execution is shown in Fig.23. 
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In Fig. 23, by uniting the battery and the motor, the motor voltage 
V M and the current I M form the observation quantities. Moreover, since 
the input torque T M of the motor becomes the load torque generated 
inside the model, the input torque T M becomes the additional load T M , 
and the output angular velocity co u becomes the observation quantity. 
It is to be noted that in the idle test of the motor, the additional load T M 
can be made 0. 

The government equation for execution of the function model 
therefor shown in Fig. 23 is given as follows: 



In equation (35), the lst-2nd rows form the state equations, and 
the 3rd-5th rows form the output equations. It is to be noted that the 
equivalent admittance Y 0 and the current source I 0 are given by the 
above-mentioned Eqs.(6) and (7). 

Moreover, the process by which Eq.(35) is derived from the 
above-mentioned Eq.(5) is as follows: 

Firstly, the 1st and 2nd columns including the moment of inertia 
J M and the inductance Lm are eliminated, and the dependent variables 
at the 1st and the 2nd rows are made x M and x L Then, the input 
torque T M of the independent variable in Eq.(5) can be substituted for 
the additional load at the 1st row the 3rd column in Eq.(5). 
3 1. 1. 2 Inertia rotation test model of motor 

With the motor, the motor voltage V M , the current I M , the angular 
velocity co M of the rotation system, and the torque T M can be measured. 
Among them, the electric system particularly has a characteristic that 
data with high accuracy are easily measured. However in the 
functional model shown in Fig.23, the motor voltage V M nearly equals 




M 



•Eq.(35) 
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the electromotive force E 0 of the battery, so that it can not be used as 
the data for identification. 

However, as for the induced voltage V M generated when the motor 
power supply is turned off, the test data for identification can be 
measured by the inertial rotation test. The model of this inertial 
rotation test is shown in Fig. 24. 

Fig.24 shows a model of the inertial rotation test from the 
moment when the power supply of the motor in the idle state is turned 
off until the motor is stopped. It is to be noted that the insulation 
resistance R c is omitted because the current is minute. The motor 
model of the inertial rotation test shown in Fig.24 is given as follows: 
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In equation (36), the lst-2nd rows form the state equations, and 
the 3rd-5th rows form the output equations. It is to be noted that the 
transient test model of Fig.23 is used for the example of the transient 
internal characteristic value identification performed here, instead of 
the inertial rotation test model. 

3. 1. 1. 3 Actual load t ransient test model of motor having operation 
system 

(1) Operation switch 

The start, the stop, the regular rotation, the reversed rotation of 
the motor, and the basic operation of the electric brake can be 
performed by a switch circuit shown in Fig.25. The test model in which 
the modeled switch circuit is incorporated into the actual load 
transient test model of Fig.21 will now be examined. 

The motor operation model can be incorporated between the 
voltage V M and the current I M which connect the battery and the motor 
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in the motor actual load transient test model shown in Fig. 21. By 
incorporating the operation switch, the actual load transient test 
model in Fig. 21 can reproduce the start, the stop, the regular rotation, 
the reversed rotation of the motor, and the drive state of the electric 
5 brake by the switch operation. It is to be noted that the simulation 
results of the model will be described after the identification of the 
transient state. 

Fig.25 is a circuit diagram of the operation switch incorporated 
between the battery and the motor. The functional model of the switch 
10 circuit is shown in the next Fig.26. 
F The resistance R L in Fig.26 is the braking resistance for the 

5 electric brake. Moreover, neither the contact resistance nor the 

m insulation resistance of the switch are considered. 

^ The operation switch in Fig.26 is composed of S ws for starting and 

=F 15 stopping the motor, S WL for operating the regular/reversed rotation, 
U and Swb for selecting necessity/unnecessity of the electric brake. As for 

K an operation signal for operating the switches, ON/OFF of the switches 

jz are respectively represented by "1" and "0". The functions of the 

^ switches are as follows: 

20 ® The switch S ws supplies the power to the motor when S ws =l, and 
stops the power supply when S ws =0. 

(g) As for Swl which controls the regular/reversed rotation, a 

switchover switch variable Swl q assumes "1" when S^=0, so that the 

current I M of the electric circuit in Fig.23 flows in the direction of the 
25 arrow to rotate the motor regularly. An S WL1 assumes 11 1" when Sw L =l, 

so that the current flows reversely to rotate the motor reversely. Also, 
o and Sw L _i are the switchover switch variables prohibiting that 

both of S WL _ 0 and Sy^.x simultaneously assume the state of "1" or "2". 

On the functional model, this operation makes S WL on the current side 
30 switch over the polarity of the output current I s , so that S WL on the 

voltage side restores the input voltage V s which has come back to the 
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original polarity. 

(3) As for the braking by the electric brake, the electric brake is made 
ON when Swb-1 and OFF when Swb=0 under the condition of the 
power supply being made OFF, i.e. S ws =0. The electric brake consumes 
5 the induced current caused by the motor idling when the power supply 
is made OFF by the braking resistance R L to generate the braking 
torque. This braking torque is determined by the braking resistance R Ll 
so that the maximum braking torque is generated when the resistance 
value is 0 

10 The government equation determined from Fig. 26 is given as 

follows: 

PmI ^-S m (l-S wS ) \V M ~\ --Eq.(37) 

Si 

Ul Swe of equation (37) is a power supply switch variable for 

T switching over the motor start/stop, and regular/reverse direction of 

nj 15 the current flow, which is expressed by the following equation: 
™ S WE =S fF5 {s WL0 -S VLl ) -Eq.(38) 

Q The government equation of the power supply system obtained by 

mutually substituting the voltage V B and the current I B of the 
government equations of the battery and the operation switch in 
20 Eqs.(4) and (37) to be integrated is given as follows: 



fa 



Sw A in equation (39) is assumed to be a braking switch variable by 
which ON/OFF of the power supply and a necessity/unnecessity of the 
electric brake are selected, and is expressed by the following equation: 
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^wvi ~ ^ ws + ~TT~S H s B (1 ) 



■Eq.(40) 



In equation (40), as for Swe 2 in the upper equation, the switchover 
switch variable member of the regular/reversed rotation assumes 
(Swl o-Swl_i) 2 =1, and S ws 2 assumes either "0" or "1 M from equation (38). 
Therefore, the upper equation can be expressed by the lower equation 

arranged by Swe 2 =S ws . 

Finally, the government equation obtained by mutually 
substituting the motor model of the actual load transient test model in 
Eq.(32) and the voltage V M and the current I M of and the power supply 
system in Eq.(39) to be integrated is expressed by the following 
equation: 



20 



E • 


"0 
0 






0 


0 




M M 














0 




0 






0 








ru 
ru 
i== 


0 

&l 




0 
0 


0 
0 


1 

O 


*G 

Kg 


0 
0 


0 
0 


0 
0 


+ 


x l 

*c 

X M 


s_ s 






0 


0 


0 


1 


0 


0 


0 




X L 




v? 




0 


0 


0 


0 




0 






X C 








0 


0 


0 


0 


Rc Y 0 S fTA 


0 


0 




1 




Is . 




0 


0 


0 


0 


0 


N G 







— Eq.(41) 

In equation (41), the lst-3rd rows form the state equations, and 
the 4th and after rows form the output equations. The abstracted 
equivalent admittance Y 0 and the current source I 0 in Eq.(41) are 
expressed by the following equations: 

1 -Eq.(42) 



R B + R c S m 

-Eq.(43) 

As an applied example of operating the motor by the electronic 
circuit incorporated into a control unit or the like, it is possible to 
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incorporate the electronic circuit into the actual load transient test 
model in Fig. 21 in the same way as the operation switch circuit of the 
motor in Fig.25. 

Also, it is possible to replace the operation system of the motor 
including the above-mentioned switch mechanism by the block 
replacement method. It is to be noted that while there is another 
method by which the electronic circuit is incorporated by using high 
circuit technology, the minimum function necessary for driving the 
motor is modeled here. Accordingly, the non-linear characteristic of a 
semiconductor, the transient response of the electronic circuit and the 
like are excluded. 
(2) Switching method 

By the switching method, the ON/OFF ratio by which the electric 
power is supplied at a cycle which does not influence the motor speed is 
changed to perform a control. The switching method can be realized by 
periodically switching over S ws in Fig. 26. The electronic circuit 
diagram close to this method is shown in Fig. 27 and its functional 
model is shown in Fig.28. Figs. 27 and 28 respectively show a model 
including the battery. 

In Fig.27, a signal voltage E CNT of a rectangular wave applied 
from the outside is converted into a base current by a resistance R A of a 
transistor T R to make the transistor T R perform an ON/OFF operation. 
At this time, E 0 and R B respectively denote a battery voltage and an 
internal resistance. D x in Fig.27 denotes a wheeling diode, and has the 
following two functions: One is a function of restoring the energy saved 
by the inductance of the motor winding generated when the transistor 
is turned OFF to the motor, and another is a function of guarding the 
transistor from the high voltage induced by the inductance. It is to be 
noted that a mechanism model of a switching signal generator is 
omitted. 

In the switching operation of Fig.28, a conductive state of ON is 
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represented by a low resistance Ron and an interrupted state of OFF is 
represented by a high resistance R 0 ff> in which both are switched over. 
Also, the function of the wheeling diode is represented by a condition 
judgment Swd shown by a hexagon frame and a resistance R D of the 
conductive state. Namely, the condition judgment S WE judges the 
voltage E CNT of the signal generator with the detected voltage E TH , the 
conductive resistance R Q n is selected when Sw E =l> and the interception 
resistance R 0 ff is selected when Selection is performed by the 

switchover switch variables S V ve_i an ^ Swe_o- 

Moreover, as for the accumulated energy of the inductance, the 
release period is judged by the counter electromotive force of the 
inductance which has turned to be negative, so that the energy is 
restored to the motor through the conductive resistance R D of the 
wheeling diode. The functions of the braking resistance R L shown in 
Fig.26 and the resistance R D in Fig. 28 appear to be the same at a 
glance. However, the former functions by a positive induced voltage of 
the motor, while the latter functions by a negative counter 
electromotive force by the inductance. The government equation of the 
switching method taking these points into account is given as follows: 



It is to be noted that the switching signal generator is excluded. 
(3) Fixed current control method 

The electronic circuit for controlling the motor by continuously 
changing the motor current or voltage is as follows: 

The basic function of the motor is to convert the given current 
into the drive torque by the torque coefficient % T . Accordingly, when an 
electronic control is performed to the motor, it is a basic function of the 
fixed current control to freely set a current targeted and to maintain 
the current. Fig. 29 shows one of control circuits for maintaining the 
output current I s in proportion to the value of a designated current 




In Fig. 29, a triangle shows an operational amplifier or the like 
which has a high amplification ratio inside. The current I CNT targeted 
is converted into the voltage by the resistance R s to be inputted to + 
5 side of the amplifier. Also, the controlled current I M is converted into 
the voltage by the current detecting resistance R E to be returned to - 
side of the amplifier, thereby performing a feedback control. It is to be 
noted that E 0 and R B are respectively the induced voltage and the 
internal resistance of the battery. 
10 Additional brief examination of the fixed current control will be 

described as follows: 
2 The voltage applied to + and - sides of the amplifier are 

r± respectively Vajnp^RsIojt and V A ^ IN p2=R E I M , as seen from Fig.29. Also, 

assuming the amplification ratio of the amplifier is A^, the output 
%1 15 voltage Va out is given as follows: 

V A OUT = (^S^CNT ~ ^E^m)^MP tti / a r*\ 

—Eq.(45) 

Also, the amplifier in Fig.29 has a high amplification ratio. 
Therefore, if equation (45) is arranged assuming 1/A MP =0, the equation 
of the output current I M is given as follows, where R S /R E in the 
20 equation indicates the current amplification ratio: 

Im^Iott -Eq.(46) 

V a _out of Eq.(45) cannot generate the voltage equal to or more 
than a supply voltage V B shown in Fig.29. Also, a voltage difference of 
an offset voltage V G fs is necessary between an output voltage upper 
25 limit value of the amplifier and a power supply voltage. This 
relationship is expressed by the following inequality: 

Va_out$V l -V ofs -Eq.(47) 

The output voltage V A _ OUT of the inequality (47) restricts the 
upper limit value of the output current I M , indicating the limit of the 
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fixed current control. The output current I S l exceeding this limit is 
determined by the closed circuit of the voltage of the battery induced 
voltage E 0 and the offset voltage E OFS shown in Fig. 29, as well as the 
voltage drop by the motor winding resistance R M , the current detecting 
resistance R E , and the battery internal resistance R B . However, since 
E 0 is not included in the functional element of the fixed current control, 
it is necessary to replace E 0 with the supply voltage V B which can be 
added to the functional element. The current I S l considering this point 
is expressed by the following equation: 

7 « s g"fr -Eq.<48) 

The current consumed by the internal circuit of the amplifier is 
very low. Therefore, by neglecting the current, equation (48) turns into 
the approximate equation of the current I B supplied by the battery. 

From the above-mentioned description, the fixed current control 
circuit in Fig.29 has a function of switching over the controlled current 
Im given in equation (46) to the uncontrolled current I SL given in 
Eq.(48) to be supplied. The functional model of the fixed current 
control in which the above-mentioned contents are arranged and 
modeled can be shown in Fig.30. 

Firstly, the government equation of the fixed current control can 
be obtained from Fig.30 as follows: 



Yb. 



0 Swamp _ \ ~7T~ ^wamp_o^cnt 

E 

1 i^E°WAMP_\ D WAMP 0 CAT 



•Eq.(49) 



In equation (49), S WAM p_i and S WAMP _ 0 indicate switchover 
elements with the condition judgment for judging the limit of the fixed 
current control. 

Finally, as for the fixed current control provides an in-control 
while the current I B within the control range of Eq.(46) is equal to or 
less than the current I SL exceeding the control limit of Eq.(48), and a 
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non-control when the current I B exceeds the current I SL . This condition 
judgment Swamp provides the following equation. It is to be noted that 
/ B and I S L in equation (50) are estimated state values. 

ifO B <U) &™S WAMP = l elkse S WAMP =0 ••■Eq.(50) 
5 In Eq.(50), S WAMP =1 means that the condition is fulfilled so that 

the current is directly supplied from the battery, while S wa mp~0 means 
that the condition is not fulfilled so that the fixed current control is 
performed. 

3. 2 Determination of initial characteristic value 
10 Characteristics which have not yet been clarified are mainly 

p transient characteristic values inside the motor. In order to 

P understand the relationship between the motor function and the 

SS functional model, the characteristic value is initialized by providing 

ifi the design value shown in the following Table 4 or an initial 

j 5 15 characteristic value obtained from a similar part, by applying a 
steplike voltage to the motor in the idle state, and by sequentially 
updating the characteristic value experimentally It is a matter of 
course that the process and method performed here can be simulated 
by the computer. 

TABLE 4 



CHARACTERISTIC NAME 


SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


INERTIA MOMENT 




[Nm/(rad/scc 2 )] 


0. 4X10' 6 


INDUCTANCE 


**» 


[H] 


3.0x10-3 



20 

Upon the identification, the simulation results at the time when 
the functional model, corresponding to model when T M =0 in Fig. 2 3, in 
the idle state is initialized by providing the characteristic shown in the 
Table 4 to Eq.(35) are shown in Figs.31(l) and 31(2). It is to be noted 
25 that for the steady internal characteristic value and the test condition 
except those shown in the Table 4, the above-mentioned identification 
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results of the steady state shown in the Table 3 are used. Also, since 

the additional load T M is idle, it is made "0". 

Fig. 3 1(1) shows the motor current I M [A] and the angular velocity 

co M [rad/sec]. In Fig.31(l), a thin dotted line indicates natural data of a 
5 measured current, a thick dotted line indicates a current I m _meg [A] 

where signal processing is performed to a measured result described 

later with a digital filter, a thick solid line indicates a current Imjsim [A] 

of the simulation result, and a thin solid line indicates an angular 

velocity co M [rad/sec] of the motor. 
10 Also, in Fig.31(2), a thick solid line indicates a current deviation 

AI ERli [A] of an initial value obtained by subtracting the current I m ,meg 
O after being signal-processed with the filter from the current I m _sim 

determined by the simulation, and a thin dotted line indicates a 

current deviation AI P err [A] at the time when a moment of inertia J M 
%1 15 described later is identified. 

s In Fig.31(2), the current I M shows a rise immediately after the 

fU power supply is turned ON where a high starting current flows, and 

Li decreases with a lapse of time, being stabilized to the idle current. 

Oppositely, the angular velocity co M rises with a lapse of time to 
20 stabilize in the idle rotation number. The change of the motor current 
indicates that the induced voltage of the motor is governed by the 
angular velocity co M through the motor constant M M . 
3. 2. 1 Noise process of measured current 

The motor current I M is measured by connecting a distributary 
25 resistance R s of 50 [mQ] between the motor and the power supply. The 
starting current measured at this time is shown by a thin dotted line 
which varies extensively in Fig.31(l). The extensive variance indicated 
by the actually measured value is due to the brush and a commutator 
of the motor, and is remarkably shown in a low rotation region 
30 immediately after the start. Particularly, the extensive variance which 
occurs at the rise of the rotation depends on a rotation angle of the 
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commutator for the brush, which is remarkably shown in a small 
motor having a few contact pieces of commutator. 

The result of performing the signal processing to the measured 
data with the digital filter in order to remove the extensive current 
5 variance is shown by the thick dotted line in Fig.31(l), for the filter 
characteristic, a low-pass filter with a cut-off frequency of 100 [Hz] is 
used. It is to be noted that since the phase delay of about 72 [msec] 
from the measured data has occurred due to filtering, the result of the 
signal processing in Fig. 3 1(1) is made coincident with the measured 

10 data by forwarding the rising point by the phase difference. Thus, the 
diagonally rising component having moved before 0 [sec] shown in 
Fig.31(l) is neglected. 

The identification processing is performed by using the data to 
which the filter processing is performed. However, it is expected that 

15 an error greatly influences the identification result since the current 
variance occurring at the start of the motor rotation is not completely 
excluded. Since this error is an uncertain variance dependent on a 
positional relationship between the brush and the commutator at the 
starting time described before, the identification is performed from this 

20 view point. 

3. 2. 2 Time history sensitivity of characteristic (process 38) 

It has been clarified that all of the characteristic values which 
form the subject of the identification are included in the motor current 
I M obtained from the output state value of the output equation in the 

25 5th row of Eq.(36). Therefore, in order to perform the identification, it 
is necessary to form a hypothesis by fully grasping the relationship 
between the characteristic values which form the subject and the 
motor current, and to verify it through the interpretation of the 
identification result. Accordingly, the folio wings are hypothesized for 

30 the influence of the characteristics on the motor. 

® It is supposed that the moment of inertia J M is proportional to the 
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change of the angular velocity (angle acceleration), provides strong 
influences during the acceleration when the change of the angular 
velocity is large, and does not provide influences in the steady state 
when the change is little. 
5 (2) It is supposed that the inductance Lm is proportional to the change 
of the motor current, provides influences when the current variance is 
large at the starting moment, and does not provide influences in the 
steady state when the current becomes stable. 

(3) It is supposed that the viscous resistance coefficient 5 M is 
10 proportional to the angular velocity, provides influences as the load of 
the internal loss with the increase of the angular velocity, and assumes 
= a fixed load in the steady state. 



H 1 ® It is supposed that the kinetic friction torque F MF is a constant and 

so is a characteristic not influenced by the state value change. The 



Ul 



15 motor has a structure of reducing the friction torque as much as 
L possible, so that the influence is little. 

i 

It is seen from the above-mentioned contents that the respective 
ways in which the characteristic values contribute are different in the 
process of a time transition from the starting moment to the steady 
20 state. Accordingly, if the sensitivity of the characteristic values for the 
change of the motor current I M is supposed to be the index of the 
identification, it can be easily presumed that the identification is 
possible. It is to be noted that the above-mentioned (3) and (D are 
added for examining the time history sensitivity of the steady internal 
25 characteristic value. 

The sensitivity of each characteristic value will be examined 
based on the above-mentioned result as follows (see flow chart in 
Fig.32): 

It has been described that the simulation result in Fig.31(l) has 
30 been performed at the initial values of the characteristic values, and 
the current deviation AIerr at that time is shown in Fig.31(2). In 
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order to make this the index of the identification, the influence which 
the characteristic value variance provides to the current deviation A 
I ERR is clarified, so that the relationship between the characteristic 
value and the motor current can be grasped. This current deviation A 
5 Ierr can be expressed by the following equation (at step S3), with a 
current measured value Imj^eg (at step Si in Fig. 33) and a simulation 
value I m _sim ( a t step S2) being made a deviation: 

A^ERR = Im_SIM ~ I\f_MEG ' ' *Eq.(51) 

If the deviation A 1^ lies within an allowable range, the 
10 transient identification is completed after verifying the identification 

H result by the simulation described later (at steps S4 and S10). 

Q 

O From equation (51), the value of the simulation is high with 

m respect to the measured value when the current deviation AI^ of 

jj? Fig. 3 1(2) is positive, and otherwise low. Moreover, an extensive violent 

=f- 15 variance is included at the starting moment. Accordingly, in order to 

H* use the current deviation AI ERR as a sensitivity of the motor current I M 

pj for the characteristic change, it is required to remove the variance, 

g Also, the offset quantity by the current deviation AI ERR is included in 

^ the sensitivity. 

20 A variance current deviation A I P err of the characteristic 

representing the sensitivity is expressed by the following equation to 

remove the variance portion of the deviation: 

err-^P SIM ~ ^ M MEG ~ ^ ERR\ _ 

" " " \ --Eq.(52) 

— 2 P SIM M SIM ) 

On the upper term of equation (52), the I P _ sm (at step S5) 
25 indicates the simulation value at the time when the transient internal 
characteristic value is changed, and the initial value current deviation 
A Ierr indicates the current deviation of the measured value Im_meg for 
the simulation value Imjsim of the initial value by Eq.(51). 

From this equation, the variance current deviation AI P ERR of the 
30 characteristic value can use the deviation of the simulation value I P _ SIM 
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of each characteristic value variance for the simulation value I M SIM of 
the initial value as a sensitivity as shown in the lower term of Eq.(52). 
However, since an identification accuracy is exceedingly evaluated by 
the initial value current deviation AIerr, it is necessary to proceed the 
5 identification while the AI^ is contained within an appropriate range 
at any time during the identification. 

3. 2. 3 Identification method by time history sensitivity change 

Firstly, the result (at step S6) upon having determined the 
current deviation AI P ERR at the time when the characteristic values 

10 are changed by +10 [%] and — 10[%] with the time history is shown in 
Figs. 33B-33E. The four diagrams show the sensitivity seen with the 
time history at the time when the characteristic values are changed. 
Fig. 33(1) is a diagram for comparing the sensitivities at the time 
indicated by thin dotted times in the lower four diagrams, so that the 

15 time history data are stored. 

Also, in the respective current deviations of Fig.33(2)-33(5), J 
indicates the moment of inertia, L the inductance, D the viscous 
resistance coefficient, and F the kinetic friction torque, which are 
generally referred to f, P (parameter)" as mentioned above. However, 

20 since the vertical axis scale in each diagram is the same as that in each 
diagram on the lower term, it is impossible to compare the four 
diagrams in parallel. Also, the solid line shows a result upon having 
increased the characteristic of the Table 4 by +10 [%], and the dotted 
line shows a result upon having decreased the same by — 10 [%]. This 

25 calculation is performed to all of the characteristic values of the 
identification subject (at step S7). 

It is seen that the characteristic becomes linear, since the result 
upon having changed the characteristic value by ± 10 [%] is 
symmetrical with respect to positive and negative in Fig.33(l). The 

30 moment of inertia J M and the viscous resistance coefficient 6 M become 
directly proportional to the motor current, and the inductance Lm 
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becomes proportional to the motor current in which the positive and 
negative are inverted. These relationships almost agree to the above- 
mentioned hypothesis of the relationship between the motor current I M 
and the characteristic, so that the variance of the deviation is removed. 
5 The identification is performed with the current deviation A 

Ip_err at the time when the characteristic of Figs.33(l)-33(5) is changed 
being made a time history sensitivity. In the identification, trials are 
repeated so that the current deviation lies within an allowable range 
by the following process (at step S4): 
10 Firstly, since the positive current deviation is large in the range 

of 0.02-0.06 [sec] from the initial value current deviation AIekr of 
o Fig.31(2), the moment of inertia J M in Fig.33(2) whose sensitivity is 

58 high in this region as shown in Fig. 33(1) is selected (at step S8) to be 
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reduced (at step S9) in the identification of the characteristic values. 
15 Since a variance value DI ERR at the starting moment comes near the 

H 3 upper side as shown in Fig. 3 1(2), the inductance Lm of Fig. 33(4) is 

nJ 

HI increased. 

g As for the order of identifying these characteristic values, the 

above-mentioned steps S2-S9 are sequentially executed from the value 

20 having larger coefficient with the mutual relationship between the 
initial value current deviation A and each variance current 
deviation A I PJERR being made an index. Alternatively, when it is 
thought that the variances do not interfere with each other due to the 
time difference as shown in Figs. 33(2) and 33(4), the identifications 

25 may be simultaneously performed by the steps S2-S9. 

It is to be noted that since the characteristic curve of the moment 
of inertia J M is similar to the current deviation AI ERR shown in 
Fig.31(2) in the curve of the time history change, only the moment of 
inertia J M may be changed for the identification under the judgment 

30 that the current deviation A Ierr is most influenced by the moment of 
inertia J M . 
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Among above-mentioned processes, the process firstly performed 
in which the updated value of the moment of inertia J M is determined 
is shown in Figs.34(l) and 34(2). 

Fig.34(l) is a diagram showing the current deviation AI ERR in 
5 Fig.31(2) up to 0.06 [sec]. Point A in Fig.34(l) indicates a deviation 
quantity a which is desired to be reduced and the position (point of 
0.03 [sec]). Fig.34(2) shows a sensitivity of the moment of inertia J M at 
the point of 0.03 [sec] shown in Fig.33(l), in which the sensitivity point 
0 (abscissa) is adjusted to the point A in Fig.34(l). 
10 It is to be noted that the current deviation shown by the both 

ordinates is adjusted to the same scale. As for the deviation quantity a 
=0.2 [A] which is desired to be removed from Fig.34(l), it is seen from 
Fig. 34(2) that the sensitivity becomes negative. 

Accordingly, a correction ratio j3 =16.8 [%] of moment of inertia 



.- ^= 
HI 

HF 15 J M is determined from an intersection B between the abscissa of the 

M= deviation 0 in Fig. 34(1) and the extension of the negative sensitivity in 

m 

nj Fig.34(2). From the Table 4, the updated moment of inertia J M assumes 

0.4 X 10* 6 X (1-0. 168) =0.3328 X 10 6 . 

The result upon having updated the moment of inertia J M is a 
20 thin dotted line showing the variance current deviation AI PJERR of 
Fig.31(2). The above-mentioned process is expressed by the following 
equations (at step S9): 
cc(t) 

♦Eq.(53) 



O 



M0 



In equations (53), the upper term is an equation of the correction 
25 ratio j3 of the characteristic value, and the lower term is an equation 
of the characteristic value after the update. As for the variables in 
Eqs.(53),a(t) [A] is sl current deviation quantity of the initial value 
current deviation A I ERR at a time "t", l(t) [A/%] is a sensitivity of a 
characteristic value at a time "t" similarly, P m _new is a characteristic 
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value after the update, and Pm_old is a characteristic value before the 
update. 

When the motor current and the characteristic value have a non- 
linear characteristic (time history sensitivity becomes asymmetrical) 
which are not proportional to each other, and the change range of the 
characteristic value greatly exceeds ± 10 [%], it is necessary to again 
determine the sensitivity from Eq.(52) to repeat the identification. 
3. 3 Result of transient identification and simulation (at step S10) 
(1) Identification result and simulation thereof 

The characteristic values via the above-mentioned process of the 
transient identification are shown in Table 5. As test data for 
identification, a result of an idle test is used in which the setting 
torque of the fixed torque transient test is assumed to be 0. 

TABLE 5 



CHARACTERISTIC NAME 


SYMBOL 


UNIT CHARACTERISTIC VALUE 


INERTIA MOMENT 
INDUCTANCE 


J* 


[Nm/(rad/sec 2 )] 0. 335 x 1 0" 6 
[H] 4.2x10-3 



The simulation result in which the transient identification result 
is reflected is shown in Fig. 35. The simulation has been performed by 
providing the transient identification result of the Table 5 and the 
characteristic value of the steady identification result of the Table 3 to 
the government Eq.(35) for execution of the fixed torque transient test 
model shown in Fig. 23. 

In Fig.35(l), the result out of consideration of the inductance Lm 
is added by a thin dotted line and the natural data of the motor current 
are omitted, while the others are the same as Fig.31. It is to be noted 
that the simulation out of consideration of the inductance has been 
performed by the same Eq.(35) with Lm-10" 12 [H] being made a minute 
value without influence. 

When both are compared, it is seen that by the influence of the 
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inductance Lm the starting current which steeply rises at the starting 
moment is saved as an energy to be released with the decrease of the 
current. 

As a result, the rising of the motor current and the peak value are 
5 suppressed in a solid line (Lm>0) with respect to the thin dotted line 
(Lm 0), and the falling part shifts to the right side. Also, the 
identification results are almost coincident as shown by the current 
deviation A 1 ERR of Fig. 35(2), while the variance at the starting moment 
appears as it is in the form of the current deviation as mentioned 
10 before. 

(2) Simulation of actual load transient test model of motor 
□ Since this identified motor model becomes a virtual prototype 

h which accurately reproduces a motor unit, a model of a virtual test 

^ which reproduces the driving state of the actual machine can be made 

=C 15 by incorporating the functional model for operation in Fig.26 into the 
M= actual load transient test model shown in Fig.21 (see Fig.42 described 

nj later). The simulation result will be described: 

5? The simulation is performed by providing the characteristic value 

M= of the electrically driven arm of the following Table 6 to Eq.(41). 

TABLE 6 



CHARACTERISTIC NAME SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


DECELERATION RATIO N Q 




1/200 


OPERATING ANGLE 9 S 


frad] 


1.57 


STOPPER STIFFNESS K s 


[N/m] 


7.0X10 4 


STOPPER RADIUS L s 


[m] 


0. 04 



20 

It is to be noted that the characteristic values other than those in 
the Table 6 are the same as those of the motor model after the 
identification. Also, the motor operation is performed by the power 
supply switch S ws , the regular/reversed rotation switch Swl> ar *d the 
25 braking switch S W b shown in Fig.26. 
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Firstly, the simulation result of the motor start and stop is given 
as follows: 

The power supply switch S ws is turned ON to start the motor 
with the regular/reversed rotation switch S WL being made OFF and the 
5 braking switch Swb being made OFF. Then, the load of 2 [N-m] is 
provided after 120 [msec], and the power supply switch S ws is turned 
OFF after 280 [msec] with the load being provided. The simulation 
result is shown in Figs.36(l)-36(4). 

Figs.36(l)-36(4) sequentially show a motor current I M [A], a load 
10 torque T L [Nm], a motor voltage V M [V], a motor angular velocity co M 
m, [rad/sec], and an operation angle of arm 0 R [deg]. The motor current 

~ I M increases, and the angular velocity co M decreases from the moment 

p the load is added. At this time, the change of the operation angle 0 R 

S 1 becomes rather slow. 

yi 

^ 15 Also, when the power supply is OFF, the motor current I M 

L assumes 0 [A], the angular velocity cu M decreases, and the load drives 

J}{ the motor from about 310 [msec], so that the reversed rotation occurs. 

H; Also, as for the motor voltage V M after the power supply is turned OFF, 

U the inductance Lm generates a spike-like voltage at the moment the 

20 power supply is OFF, and then the induced voltage V w by the angular 
velocity co M of the motor occurs, so that the state where the motor is 
changed into the generator is described. 

The simulation result of the motor regular/reversed rotation, as 
well as the braking is as follows: 
25 Figs.37(l)-37(4) show an example in which the motor is stopped 

with and without the electric brake by OFF/ON of the braking switch 
S VV B- Figs.37(l)-37(4) show states as follows: The motor is started with 
the regular/reversed rotation switch S WL being made OFF and the 
power supply switch S ws being made ON to be compulsorily stopped 
30 with the stopper. Then, the motor is rotated reversely in 0.6 [sec] after 
the regular/reversed switch Swl being made ON and is again 
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compulsorily stopped with the stopper on the opposite side. After the 
lapse of 1.4 [sec], the power supply switch S ws is made OFF, so that the 
motor is stopped. 

Figs.37(l)-37(4) sequentially show a motor current I M [A], an 
5 operation angle of arm 6 R [deg], a motor angular velocity co M [r ad/sec], 
and a motor torque T M [Nm]. Also, the solid line of Figs.37(l)-37(4) 
shows the state where the electric brake by the braking resistance of 
R L =1 [m£2] is applied, while the dotted line shows the state without the 
brake. Both of the regular and reversed rotations generate an 
10 attenuated oscillation at the time of collision at the stopper. 
H- Also, when the motor is compulsorily stopped by the stopper, a 

□ flexible structure of a stopper gum absorbs the maximum torque of the 

motor to be compressed. If the power supply is cut off in this state, the 

^ arm is pushed back by the compression power owing to the bend of the 

U] 

=F 15 stopper. The effect of the electric brake appears at the stop state after 

3 

M 1.4 [sec]. It is seen from the motor angular velocity co M and the 

operation angle of arm 0 R that the motor idle by the stopper 
^ counterforce is well absorbed as shown by the solid line. 

H 4. Motor identification apparatus 

20 Hereinafter, an apparatus required for identifying the above- 

mentioned test model of the steady and transient states of the 
functional parts will be described. 

The functional parts incorporated into the products have a 
structure working with assembling the functional parts to the 
25 input/output portion and by mutually associating them. As shown in 
Fig.21, the battery for driving, the electrically dx*iven arm forming a 
load and the like are necessary e.g. for the motor test. Also, in order to 
control the motor, the operation system shown in Figs. 25, 27, and 29 is 
necessary. 

30 The test model for thus performing an identification by using the 

test data of the actual machine requires modeling the functional parts 
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which form the driving source and the load of the tested object. Subject 
to the influence by other functional parts in this way, it is possible to 
standardize the function and the behavior as the specification of the 
tested object and the test standard. If there is something which 
5 operates and reproduces the behavior just as standardized function 
and behavior, it can be used as a substitute common to each test 
method. 

If the functional model which forms the substitute of the driving 
source and the load is modeled and incorporated into the test model as 

10 an identified model, the identification defined to the internal 
characteristic value of the tested object can be performed. Also, as for 
the testing unit and the test model, the group of the tested objects 
whose parts specifications are different can be replaced with block 
replacement means. 

15 Fig. 38 shows an example in which the concepts are applied to the 

motor identification apparatus. In the identification apparatus of 
Fig. 38, a motor model 1 has a block replacement structure for a test 
reproduction model 22. A reproduction model 2 is composed of a drive 
model 21 such as a battery and an operation model, a load model 22 

20 such as a load generator model, an observation model 23 for observing 
the state value of the motor model 1, and a characteristic update model 
24 for updating the internal characteristic value. 

Also, a drive circuit unit 4 takes charge of the battery and the 
operation system which forms driving source of a tested motor 3. A load 

25 generator 5 and an electric load 6 take charge of the mechanical load of 
the motor 3. The driving source and the mechanical load are modeled 
as a drive model 21 and a load model 22 already identified. 

Furthermore, as for the relationship between the actual machine 
testing unit and the models, the former drive model 21 has a structure 

30 of controlling the drive circuit unit 4 through a drive system control 
unit 71, and the latter load model 22 has a structure of controlling the 
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load generator 5 through a load system control unit 72. 

The measured value of the motor 3 is obtained at a measurement 
unit 8 to be provided to a calculation unit 9. The calculation unit 9 
receives the observation value of the observation model 23 and the 
5 measurement from the measurement unit 8 to update a characteristic 
update model 24. 

The motor drive model 21 shown in Fig. 38 corresponds to the 
model in which the battery on the left side of Fig. 21 and the operation 
system circuit of Figs.26, 28, and 30 are united. The load generator 
10 model 22 corresponds to the arm mechanism on the right side of Fig.21 
and the fixed torque generation model of Fig.23. 

The operation and the behavior reproduced on the former model 
are reproduced at the drive circuit unit 4 by the drive system control 

unit 71 to be provided to the motor 3 of the tested object. Similarly, the 

in 

V 15 load state provided to the motor 3 is reproduced on the load generator 
f. 5 by the load system control unit 72 controlling the electric load 6. 

yJ Accordingly, if the motor functional model 1 shown in Fig. 4(1) is 

H= connected to the reproduction model 2, the drive model 21 and the load 

rf model 22 preliminarily provided corresponding to the drive system and 

20 the load system turns into predetermined steady test models and 
control the drive circuit unit 4 and the electric load 6. 

Together with this operation, the steady test data generated by 
the motor 3 is measured by the measurement unit 8, so that the 
calculation unit 9 which has received the data provides the 
25 government equation of the steady functional model i.e. the steady 
internal characteristic value to the motor model 1. Thus, the steady 
identification is completed. 

Then, the reproduction model 2 changes the drive model 21 and 
the load model 22 into the transient test model for the transient 
30 identification. Together with this operation, the transient test is 
executed, so that the transient test data are provided to the calculation 
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unit 9 from the measurement unit 8. The calculation unit 9 compares 
the value obtained by substituting the already calculated steady 
internal characteristic value for the internal characteristic value 
inputted from the observation model 23 with the measured data from 
5 the measurement unit 8. The internal characteristic value of the motor 
model 1 is properly updated as mentioned above, by the characteristic 
update model 24 so that an error between them becomes small to lie 
within an allowable range. Thus, the transient identification ends, so 
that all of the identification works are completed. 

10 5 . Virtu al test 

The model which reproduces the function, the performance, and 
the characteristic value of the parts or the products composed of the 
m parts is identified with the test data of the actuality, so that the 

internal characteristic value and the reproductivity are verified. 
15 Therefore, it can be utilized as a virtual prototype for performing 
various evaluation tests on the computer instead of the actual machine. 
An example where the virtual prototype is applied to the product 
development will be described. 

R. 1 Relation ship between product development and modeling 
20 The development process of product and the formation process of 

model are shown by the relationship in Fig. 39. 

Fig. 39 shows a relationship between a development process of a 
project stage, a design stage, a prototype and test stage, and modeling 
and identification performed in parallel with the development. Fig.39 
25 shows a process including a concept of concurrent engineering in which 
development processes of the product are performed in parallel. 

On the other hand, the process of modeling has a flow as follows: 
An abstract product target at the project stage is firstly modeled, the 
model is sequentially detailed along with the development process, a 
30 validity of the project and design contents by the model is verified, and 
the identification is finally performed with the test data of the 
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prototype, so that the virtual test is performed executing computer 
simulation with the identification model instead of the various 
evaluation tests. 
(1) Relationship between product development and modeling 
5 If this virtual test is applied to the development of the car, for 

example, the process is given as follows: 

As for the abstract model at the project stage such as a model 
reproducing a car performance which forms the target characteristics 
of the projected car such as an acceleration time of 0-400 [m], a 
10 stationary fuel cost performance, and an acceleration performance, a 
M; possibility of achieving the target characteristics is verified by the 

O model assembled with the characteristics such as a weight, a 

£^ resistance, a power characteristic, and a gear ratio, 

ijj In the next product part formation and basic design, the model is 

=P 15 disassembled following the formation method of the basic functional 
parts such as an engine, a transmission, a body realizing the target 
ry characteristics, and modeling is performed with the internal 

S characteristic values of the functional parts being clarified, so that the 

function and the performance required for the functional parts, and the 
20 possibility of realizing the target characteristic value are verified. 

Similarly, in the detailed design, the modeling is detailed 
according to the contents of the detailed design of an engine, a clutch, a 
brake, a transmission, a control unit, and the like, whereby the validity 
of the design contents in the details such as a torque variance, a 
25 transmission shock, a vibration noise, and the control algorithm is 
verified with the model. 

The model made by the above-mentioned process is sequentially 
disassembled into the lower model from the abstract model at the 
project stage along the part formation, so that a systematic model 
30 hierarchized from the product to the part and its details is formed. 

Also, in the above-mentioned process of the project, design and 
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modeling, the development is proceeded by the design drawing of the 
actuality realizing the target characteristic, and the model verifying 
the design drawing. However, both are eventually hypotheses. 
Accordingly, it is necessary to test both by an experimental car made 
5 based on the design drawing or the like, to identify the model, and to 
verify the validity of the hypotheses. 
(2) Relationship between product development and identification 

By the above-mentioned reason, the identification shown in the 
lower side of Fig.39 is performed based on the data of the functional 
10 parts and a car test. The identification and the test necessary therefor 
are performed separately in the above-mentioned steady and transient 
O states. The model is converted in accordance with the test method as 

S3 shown in Figs.9-13 for the steady test, and in Figs.21-30 for the 

si 

in transient test, so that the identification of the test model is performed 

15 by the test data. In the identification performed here, the identification 
^ result of the test models firstly performed per functional part are 

fU integrated along the part formation to be made the identification of the 

b product model. 

It is supposed that the functional parts of the battery, the motor, 
20 and the electrically driven arm in Fig.21 are individually modeled, for 
example, and the product in which the parts are mutually connected is 
an electric actuator. The model of the electric actuator can be 
represented by a hierarchized model where the models of the battery, 
the motor, and the electrically driven arm are made lower level. Also, it 
25 is seen that the equation in which the models are integrated turns into 
Eq.(34). 

Accordingly, if the motor at the center of Fig.21 and the 
electrically driven arm on the right side of Fig.21 are individually 
identified, and they are made the characteristic value of Fig.21, the 
30 electric actuator on the upper hierarchy is naturally identified. The 
model which can faithfully reproduce the characteristic value and the 



64 



performance of the product along the mode of the actuality is called a 
"virtual prototype". 

As an applying method of this hierarchized model, it is possible to 
convert the electric actuator model into the test model when the 
5 battery and the electrically driven arm are already identified in Fig.21, 
to make the identified internal characteristic values fixed values, and 
to perform a partial identification defined to an unfixed internal 
characteristic value of the motor model. 
5. 2 Virtual test by virtual prototype 
10 Thus identified part model is confirmed as a virtual prototype of 

M* the part reproducing the function and the performance held by each 

C part. Similarly, the product model in which the part models already 

jr| identified are mutually connected and integrated means that the 

,~ model is confirmed as a virtual prototype reproducing the product. Also, 

=P 15 these models result in the identification by converting the models into 
N= test models in accordance with the above-mentioned test methods, 

fjj comparing and evaluating the reproducibility and the test result of the 

E actuality, and verifying the difference between both. 

Accordingly, the virtual prototype of the identified product and 
20 parts can be evaluated as to the validity by performing the computer 
simulation instead of the various tests of experimental articles and the 
like. Performing the evaluation test of the actuality on the computer by 
utilizing the virtual prototype is called a "virtual test" of the parts or 
the product. The simulation which has verified the validity of the 
25 identification result obtained in Figs. 37-38 is also one example. 
(1) Internal formation of virtual test model 

The virtual test is usually applied with the test standard 
determining the same environment condition as the evaluation test of 
the actuality, a test condition and performance such as driving 
30 operation, and an evaluation standard such as a success/failure 
judgment, a performance prediction, a presence/absence of abnormal 
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phenomena, and the like by the reproduction data of the performance 
and the behavior. This can be typically shown in Fig. 40. 

Fig. 40 is a basic arrangement of a virtual test model 10 where a 
virtual prototype 11 and a test standard group 12 applied thereto are 
5 modeled to be incorporated. The virtual prototype in this arrangement 
comprises two kinds of functional part models, an already identified 
functional part model and a functional part model in which design 
values and characteristic values of similar parts are provided upon the 
virtual test. 

10 Both are distinctive in that the former is used as a model of a part 

having a basic function governing the function and the performance of 
the product, while the latter is used as a model extended for various 
products whose characteristic values are different in the application, 
the structure, and the shape of the products. The product model having 

15 both models united indicates that the evaluation test of various 
products and the test evaluation of various product groups where the 
derived functional models (latter) of different characteristic values are 
combined with the identified basic functional model (former) can be 
performed by the virtual test. 

20 (2) Virtual prototype example of virtual testing car 

For example, a simple model of the virtual prototype for car in 
Fig.40 can be shown in Fig.41. Fig.41 is an example of a virtual 
prototype showing a functional model where the functional parts of the 
engine and the transmission are modeled to be integrated as a power 

25 train (P/T), and a functional model where various types of cars which 
are different in a structure and the like are modeled as a body, and 
both are integrated as a car. Between the models of the functional 
parts, the engine and the transmission are related with the angular 
velocity co c that is a potential quantity, and with the torque T e that is 

30 a flow quantity, while the transmission and the body are related with 
the torque v b that is a potential quantity, and with the torque 4 that is 
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a flow quantity. 

The characteristics in the models are as follows: 

The characteristics of the engine are a moment of inertia J e , a 
viscous resistance coefficient 6 e , and an internal torque T e . 6 e and T e 
5 among them respectively indicate a partially linearized gradient and a 
generation torque when the engine rotation number is "0", in the 
rotation — torque characteristic in which a throttle divergence a e , an 
atmospheric pressure p e , and an air temperature t e are inputted. 

The characteristics of the transmission are a gear ratio N m , a 
10 final gear ratio N d , and a stiffness of an output axis C P . Finally, the 
characteristics of the body are a tire radius R^a mass of the body M b , a 
resistance during running D b , and an ascent of a running road 0 . 
Since the characteristics of the body, Rr, M b , and D b are different 
according to the body structure or the like, they are required to be 
15 replaced according to the type of car. 

Moreover, the control characteristic values of the engine control 
unit and the transmission control unit respectively controlling the 
engine and the transmission are preliminarily identified by the test 
data. 

20 In order to perform a virtual test under the same condition as the 

car driving state by using the virtual prototype in Fig.41, the following 
model is further required. Firstly, for the gradient 0 , the atmospheric 
pressure p e , the air temperature t e and the like, the environment 
condition model provided as an environment condition for using cars is 

25 required. Furthermore, for the throttle divergence a e , and a gear ratio 
N K of the transmission, the operation model which provides the driving 
operation condition of cars is required. 

Finally, in order to evaluate the reproduction result of the virtual 
prototype in which the virtual test has been performed, the observation 

30 model for observing co e , v b , T e , f b , and the like is required. The model 
for the conditions and for performing the evaluation is modeled based 
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on the evaluation test methods of the test standard models shown in 
Fig.40. 

The virtual test performed by the model can be applied to various 
types of cars by incorporating the characteristic values of the body 
5 model which are different depending on the type of car into the virtual 
prototype of the functional part in the P/T which faithfully reproduces 
the characteristic value, performance, and behavior of the actual 
machine. The virtual test using the virtual prototype of the car model 
can be substituted for the various tests performed to various types of 
10 car on which the same P/T is mounted. 

In the development of cars so far, the evaluation tests extending 
P over a wide range have been performed to the cars having many body 

structures, so that the validity as a commodity has been verified. In 
such vast numbers of tests, the number of cars made on the 
* 15 experimental basis can be reduced and the development period 
jf_ required for the test can be shortened by the virtual test performed 

fU with the body characteristic value in the virtual prototype being 

□ replaced and with modeled test standard groups being combined. As a 

result, cost for a development can be reduced. 
20 5, 3 Virtua l testin g system (model) 

The basic arrangement of the actual machine test and the 
identification as well as the virtual test by the virtual prototype can be 
shown as in Fig. 42. 

It has been already described for the virtual test of the virtual 
25 prototype that the tested model which forms a subject of the actual 
machine test can be converted into the test model, the test standard 
can be modeled as the driving operation model and the environment 
condition model, and the execution result of the virtual prototype can 
be modeled as the observation model. 
30 It is possible to compare and evaluate both of the observation 

data reproduced by the virtual test and the test data of the actual 
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machine test and to perform the identification. This indicates that the 
virtual test model and the test model for the identification can be made 
common. 

Namely, the virtual testing system (model) can have consistency 
5 of the actual machine test, the identification, and the virtual test by 
incorporating an identification model 9 (corresponding to a calculation 
unit in Fig.38) which models the identification technique and process 
to update the internal characteristic value of the virtual prototype, and 
an evaluation model 35 in which the virtual test verifies the validity of 
10 the reproduction data. 

In Fig.42, the actual machine test data are inputted to the 

D evaluation model and to the identification model, and the result of the 

O 

H evaluation model is inputted to be related. As for both models, the 

following processes are performed: In case a problem arises in the 

^ 15 result of the success/failure judgment in the virtual test, the former 
performs the actual machine test to evaluate again as shown in a 

hi dotted line, and in case there is a problem in the identification result of 

[2 the virtual prototype, the latter can perform the identification again. 

P Main contents of each portion in the virtual testing system shown 

20 in Fig.42 are given as follows: 

® The virtual prototype 1 corresponds to the motor model in Fig.38, 
and a model which reproduces the function, the performance, and the 
characteristic value of the product and part, which is composed of a 
mechanical system model relating to the mechanism and the structure, 
25 and the control system model for controlling them. 

(2) A driving operation model 31 is a model for operating the virtual 
prototype 1 by the same driving condition as the actual machine. For 
example, there are an accelerator operation, a transmission operation, 
a braking operation of a car, and the like. 
30 (3) An environment condition model 32 is a model for reproducing the 
influence of the environment with the virtual prototype 1 being made 
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operate under the same condition as the environment for using the 
actual machine. For example, there are a running road, a temperature, 
a humidity, an atmospheric pressure, and the like when the car makes 
a market running. 

5 © An observation model 23 is a model for reproducing the data, on the 
model, at each measurement point measured in the actual machine 
test as shown in Fig. 38, so that the data which are difficult to be 
measured in an actual machine can be observed on the model. 
(5) An internal characteristic update model 24 is a model for updating 
10 the internal characteristic value of the virtual prototype based on the 
N- identification result, as shown in Fig. 38. 

p ©An actual machine test portion 33 performs a test which has been 

hj performed in the prior art, and a measurement unit 34 in the actual 

ri machine test portion 33 shows a unit for taking the test data into a 

*P 15 computer on real time basis. 

jU The above-mentioned operation can be summed up as follows: If 

Hi the driving operation condition and the environment condition are 

"K provided to the virtual prototype 1 from the models 31 and 32 

H respectively, the virtual prototype 1 provides the reproduction data 

20 obtained as a result of the simulation at this time to the evaluation 
model 35 from the observation model 23. The evaluation model 35 
evaluates, instead of the actual machine test, whether or not the 
reproduction data of the simulation result is valid (validity of 
product/part) by comparing the evaluation reference determined by the 
25 test standard with the reproduction data. 

In case of a failed evaluation, the evaluation model 35 instructs 
the virtual test portion 30 to change the driving operation condition 
and the environment condition, and to conduct the identification again 
if the deviation between the actual machine test data and the 
30 reproduction data is large. Thus, as shown in the dotted line in Fig. 42, 
the driving operation condition and the environment condition are 
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respectively provided to the actual machine test portion 33 from the 
models 31 and 32. 

Accordingly, if the deviation between the actual machine test data 
from the actual machine test portion 33 through the measurement unit 
34 and the reproduction data is large, the evaluation model 35 
instructs the identification model 9 to conduct the identification again. 
Since the identification model 9 corrects the internal characteristic 
value of the virtual prototype 1 by the internal characteristic update 
model 34 after receiving the actual machine test data, the evaluation is 
performed again at the evaluation model 35 through the observation 
model 23. As a result, the identification model 9 updates the internal 
characteristic value of the virtual prototype 1 so that the deviation 
between the actual machine test data and the reproduction data 
becomes least. 

When the evaluation ends, another modification (e.g. 
modification according to a type of car) in the virtual prototype 1 is 
evaluated. Particularly, by performing the identification to the internal 
characteristic value (control parameter) of the control model in the 
virtual prototype, the virtual prototype 1 can be tuned, according to the 
type of car, to an optimum state. 

5. 4 Actual machine test of transmission annlving virtual test 

In case the actual machine test is performed to a unit of the 
functional part of the tested object, the functional part of the drive 
system for driving the part and the functional part of the load system 
which forms the load are required. The method in which the virtual 
prototype where such a drive system functional part and a load system 
functional part are identified is incorporated into the testing unit, and 
the actual machine test of the tested object is performed will now be 
examined. 

Fig. 43 shows an example of a testing unit of a transmission 
included in the virtual prototype. 
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Fig.43 shows a unit for reproducing the drive characteristic of the 
engine and the load characteristic of the car by the drive motor and the 
load generator, and for performing the actual machine test to the 
transmission proper along the actual use state. This unit derives the 
5 performance required for the running of the car from the engine 
performance which forms the drive source, and the evaluate test about 
various performances relating to the acceleration such as an 
acceleration performance or a gearshift shock is performed by the 
actual machine. 

10 In Fig.43, the testing unit 33 is composed of the drive motor 4 for 

driving the transmission, the load generator 5 which forms the load, 
the electric load 6, and the measurement unit 8 for measuring the 
driving state. Also, to the testing unit 33 the virtual prototype 11 of the 
drive motor 4, the load generator 5 and the transmission 50 is 
43 15 connected for reproducing the driving state of the car. The testing unit 
33 comprises the drive system control unit 71 and the load system 
ru control unit 72 which faithfully reproduce the result at the drive motor 

Of 4 and the load generator 5. 

O These control units 71 and 72 provide the drive torque of the 

20 engine to the transmission 50 by the drive motor 4, and reproduce the 
load torque provided to the transmission 50 through the body by the 
load generator 5 and the electric load 6. Both models of the drive motor 
4 and the load generator 5 can reproduce the function, the performance, 
and the behavior of the engine and the body in the testing unit 33 by 
25 the identified virtual prototype 11. 

The virtual prototype in which these drive motor 4, transmission 
50 and load generator 5 are modeled is shown in Fig.44. It is to be 
noted that the operation system models 31, 38, 39, and the 
environment system model 32 inputted to the transmission model 1 in 
30 Fig.43 are the models of the test standard mentioned above. The 
former represents the model which inputs the driving operation 
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condition of the driver, and the latter inputs the environment condition 
such as an air temperature which influences the internal characteristic 
value of the transmission. Also, the observation system models 23, 37, 
and 40 are models for observing the state values inside the virtual 
5 prototypes. 

The motor model 13 of the drive system in Fig.44 is composed of 
the models of the motor, the battery, the fixed current control 
respectively in Figs. 6, 4A, and 30. It is to be noted that a part of the 
internal characteristic values of the motor model 13 are omitted for 
10 simplified description of the model. Similarly, on the load side, the load 
*f generator model 14 is composed of the braking resistance in Fig.4(l) 

O and the electric load is composed of the braking resistance on the left 

m side of Fig. 13, where a part of the internal characteristic values is 

5 S omitted. 

f 15 In Fig.44, the angular velocity co e and the input torque T e 

M= connecting the drive system motor model 13 and the transmission 

model 1 are provided as target values of the drive system control unit 
71. Similarly, a velocity v b and a drive force % of the load generator 
model 14 are provided as target values of the load system control unit 
20 72. Also, the drive torque T e of the drive motor 4 can be reproduced by 
controlling the control current I CNT with the operation system from the 
fixed current control model in Fig.30 and by adjusting the current of 
the drive motor. 

Also, the load torque of the load generator 5 can be reproduced by 
25 being changed into the braking resistance value through the operation 
system from the motor model of the braking test in Fig. 13. Also, as for 
the load torque, the fixed current control model and its unit in the 
same Fig.30 as the drive motor may be incorporated instead of the 
braking resistance. 

30 Moreover, into the testing unit 33 shown in Fig. 43, the internal 

characteristic update model 24 is incorporated in which the evaluation 
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model (see Fig. 42) provided in the calculation unit 9 evaluate a 
difference if generated between the test result of the transmission 50 
and the reproduction result of the virtual prototype 11, and the 
identification model (see Fig. 42) similarly provided in the calculation 
5 unit 9 identifies the internal characteristic value in the transmission 
model (virtual prototype) 1. 

* 5 Actual machine test ex ample of power train (PJT) applying virtual 
test 

An example will be examined in which the actual machine test of 
10 the P/T where the engine and the transmission, except the motor 4 on 
the drive side and the virtual prototype 13 and the drive system control 
unit 71 leading to the motor 4, are united is performed in the actual 
m testing unit shown in Fig.43. 

^ The example of the actual machine testing unit having the P/T 

u i 

*P 15 and the load generator is shown in Fig.45. 

U In this example, the same P/T 51 is mounted on the bodies whose 

ry structures are different shown in Fig.41, and the actual machine test is 

U performed to the various cars, so that the adaptability of the P/T 51 for 

i"* the cars having the different body characteristic values is verified. 

20 The virtual prototype 1 of the P/T 51 requires many kinds in 

which transmissions such as an automatic transmission, a manual 
transmission, and a continuously variable gear are combined in 
addition to kinds of motors such as a motor, a gasoline engine, and a 
diesel engine. Accordingly, these virtual prototypes 1 are replaced with 
25 a model replacement unit 25. 

Also, the virtual prototype 1 of the P/T 51 has a control model for 
reproducing the software of the control unit of the engine and the 
transmission, so that the virtual prototype 1 can control the model of 
the virtual prototype 1 and can directly control the actual machine of 
30 the engine and the transmission through a control input/output unit 
36. 
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Accordingly, when adaptability of the body performance 
reproduced at the actual machine P/T 51 and the load generator 5 is 
worsened, the internal characteristic value of the control model can be 
changed through an internal characteristic update unit 24. Then, the 
measure effect is firstly confirmed in the virtual prototype 1, and 
secondly the adaptability is again evaluated by the actual machine test 
to be confirmed. 

In case a low-speed region torque of engine mounted on a car 
whose body weight or running resistance is large is insufficient and the 
starting acceleration power is insufficient, the control characteristic 
value of the control model is updated through the internal 
characteristic update unit 24 in order to improve the acceleration 
performance and increase the fuel amount at starting the acceleration. 
The measure effect of the control amount change and the bad influence 
on the fuel-efficiency and the exhausted gas worried because of the 
increased amount are firstly confirmed by the virtual test in the 
simulation of the virtual prototype, and the actual machine evaluation 
test of the P/T 51 is finally performed by using the virtual testing unit. 

By this process, the tests for reproducing bad influence which 
examines the measure effect and many test items can be verified by 
the simulation of the virtual prototype in a short time, so that it 
becomes possible to define the actual test items requiring a long time 
based on the verification result, and to perform an efficient test 
operation. 

As described above, a characteristic value identification method 
and an apparatus therefor according to the present invention is 
arranged so that a functional model of a part is prepared based on a 
potential quantity and a flow quantity representing a strength and a 
quantity of energy applied to the part, a steady internal characteristic 
value of the functional model in a steady state is identified, and a 
transient internal characteristic value of the functional model in a 
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transient state is identified by using the identified steady internal 
characteristic value. Therefore, the function of all of the articles which 
govern two dimensions of the potential quantity and the flow quantity 
prescribing the energy can be modeled. 
5 Also, since the transient identification is performed to the 

internal characteristic value of the model after the steady 
identification, the steady identification without an interference of a 
transient state can be performed. Since the transient identification is 
performed based on the steady identification, there is an effect of 
10 obtaining an accurate internal characteristic value by a simplified 
identification process. Also, it is possible to faithfully reproduce 
products/parts in the steady state and the transient state, 
m Moreover, by manufacturing an apparatus which uses such a 

characteristic value identification method, it becomes possible to 
15 promptly perform the identification of the functional model of the same 
kind. 

Furthermore, the characteristic value identification apparatus 
according to the present invention is arranged such that the functional 
model having the characteristic value identified by the identification 
20 apparatus is incorporated in a virtual testing system as a virtual 
prototype, reproduction data are obtained from the virtual prototype 
by providing a driving operation condition and an environment 
condition, actual machine test data obtained by the driving operation 
condition and the environment condition are compared with the 
25 reproduction data, and a re-identification is performed depending on 
the comparison result as the occasion demands. Therefore, it becomes 
possible to omit the actual machine test performed in a development 
process of a design, prototyping, and a test, to shorten a period, and to 
reduce cost for a development. 
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